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ELECTRON WAVE INSTABILITIES IN A MAGNETOPLASMA 
CAUSED BY ELECTRON-NEUTRAL COLLISIONS* 


By John Q. Howell 
Langley Research Center 

SUMMARY 

A new class of collision-dependent electron waves is found in a non-Maxwellian 
Lorentz magnetoplasma, and it is shown that electron-neutral collisions may cause these 
waves to become unstable. The Boltzmann equation with collision integral is solved 
while assuming propagation to be either parallel or perpendicular to the magnetic field. 
Both conductivity tensors are derived and put in a form useful for numerical calculations. 
The full set of Maxwell's equations is then used to derive the dispersion relations for 
both directions of propagation. The dispersion relations are initially solved for a mono- 
energetic electron distribution function and then an isotropic distribution with a peak of 
nonzero half-width is treated. Some consideration is also given to an isotropic 
Maxwellian distribution both with and without a bump on the tail. 

An an example of propagation parallel to the magnetic field, transverse electro- 
magnetic or "whistler" waves in a nitrogen plasma are considered. A new collisional 
mode is found with real frequency near the electron cyclotron frequency, and under 
appropriate conditions it may be either convectively or absolutely unstable. Longitudinal 
or electrostatic waves propagating perpendicular to the magnetic field in a nitrogen 
plasma are also considered. A series of new collisional modes with real frequency near 
the electron cyclotron frequency and its harmonics are shown to be unstable under some 
conditions. For both directions of propagation, the electron waves can be unstable if a 
large fraction of plasma electrons have their velocity in a region where the electron- 
neutral collision frequency increases (or decreases) rapidly enough. The requirement 
on the electron-neutral collision frequency is easily met by nitrogen, mercury, and the 
Ramsauer gases. 


*The information in this paper is largely based on a dissertation entitled 
"Collisional Effects on Waves in a Magnetoplasma" submitted in partial fulfillment of the 
requirements for the degree of Doctor of Philosophy in Applied Physics, Stanford 
University, Stanford, California, June 1970. 



INTRODUCTION 




In recent years both experimental and theoretical investigations of collision- 
induced instabilities in a plasma have been undertaken. Most of this effort has been 
directed toward determining the effect of electron-neutral collisions on electron waves 
in a weakly ionized plasma where electron-electron and electron-ion collisions are neg- 
ligible. It has been found that both convective and absolute instabilities are possible if 
(1) the electron distribution function is sufficiently non-Maxwellian and if (2) the electron- 
neutral collision cross section increases (or decreases) rapidly enough with electron 
energy. The purpose of this paper is to obtain dispersion relations for waves in a weakly 
ionized magnetoplasma by treating electron- neutral collisions in the most rigorous 
manner possible, and then to solve these dispersion relations for realistic plasma con- 
ditions. There may be electron-neutral attachment in plasmas of this type but the 
resulting negative ions will not affect the waves being considered because of their large 
mass. 


Historical Review of Previous Work 

To predict the effect of electron-neutral collisions on waves in a plasma, two 
approaches have been utilized. Some workers have considered the effective radiation 
temperature of the plasma whereas others have determined the stability of waves propa- 
gating in the plasma. The latter point of view is taken in this report. 

In 1958 Twiss (ref. 1), while investigating the radiation spectra from a plasma, 
found that under some conditions a non-Maxwellian plasma could radiate more energy 
than it absorbs over a region of the spectrum. Then in 1961 Bekefi et al. (ref. 2) showed 
that electron-neutral collisions could be responsible for such enhanced emission. This 
phenomenon was investigated experimentally by Fields et al. (ref. 3) and they found a peak 
in the radiation spectra from some plasmas at the electron cyclotron frequency. Both 
their experimental and theoretical work showed that this condition occurred when the 
electron energy distribution was highly non-Maxwellian and when the electron-neutral 
collision cross section increased rapidly enough with velocity. The radiation tempera- 
ture approach was also used in references 4 to 8 where more experimental results are 
reported and in references 9 to 12 where further theoretical work is given. Durcek 
(ref. 12) used a purely quantum mechanical point of view and obtained results similar to 
those of other authors. 

Of those authors treating the collisional effects on waves propagating in a plasma, 
only those analyzing electron-neutral collisions by means of a realistic velocity-dependent 
collision model are discussed here. In 1964 Drummond et al. (ref. 13) assumed that the 
electron-neutral collision cross section was isotropic but velocity dependent and derived 
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the dispersion relation for electrostatic waves in an unmagnetized plasma. They recov- 
ered the well-known Landau dispersion relation in the collisionless limit. Derfler 
(ref. 14), as did Drummond, assumed the electrons to be colliding with infinitely heavy 
neutrals but allowed for an angular dependence in the collision cross section. In the 
limit of isotropic collisions, his dispersion relation is identical to that of Drummond. 
Derfler solved the dispersion relation for a simple case and discovered a group of poten- 
tially unstable collision-dependent modes he called "pseudo-sound waves." Suzuki 
(ref. 15), Idehara and Sugaya (ref. 16), and Shimomura and Mitani (ref. 17) investigated 
the effects of electron-neutral collisions on waves in a magnetoplasma. They all 
assumed a Lorentz gas model (that is, infinitely heavy neutrals and ions) and allowed for 
an isotropic velocity -dependent collision cross section. Suzuki considered transverse 
electric waves propagating parallel to the applied magnetic field and found a new mode 
with real frequency near the cyclotron frequency. This mode was unstable under some 
conditions and vanished in the collisionless limit. Idehara and Sugaya considered waves 
propagating perpendicular to the magnetic field and investigated both the ordinary and 
the extraordinary waves. They assumed monoenergetic electrons and found that either 
mode may be unstable at frequencies near the electron cyclotron harmonics. However, 
their results are valid only in the near infinite wavelength limit. Shimomura and Mitani 
also considered ordinary waves propagating perpendicular to the magnetic field. Unfor- 
tunately, they made an error in solving the first-order Boltzmann equation and their 
conductivity tensor is incorrect. Their solution for the first-order electron distribution 
function does not satisfy the differential equation. The error is of such nature, however, 
that their term is correct, and hence their dispersion relation for the ordinary 

wave is correct. They also assumed monoenergetic electrons and found collision- 
induced instabilities. 


Collision Models 

When solving the Boltzmann equation 


~ + v • -2L - JL(e + v x b) 
at ar m 


_af_ _ /af 
av \at 


/ coll 


( 1 ) 


for the electron distribution function f (v'jr'jt), the collision model assumed is very 
important in determining how easily the solution will be obtained. The most common 
(and the simplest) approach merely ignores collisions altogether, that is. 



( 2 ) 
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In this form the Boltzmann equation is usually called the Vlasov equation. On the other 
hand, for elastic binary collisions, the most exact collision term one can use (and the 
most difficult to work with) is the Boltzmann collision integral 


y |v - vjcr[f (v 1 ) f(v’) - f(v) F(v)]dO d 3 V (3) 

where a is the collision cross section and f(v) is the distribution function of the 
particles with which the electrons are colliding. The velocities of an electron before 
and after the collision are v and v ', and the before and after velocities of the particle 
with which the electron collided are V and V'. The differential solid angle in velocity 
space in the center of mass frame of reference is dO. For electrons in a weakly ionized 
plasma, electron-neutral collisions are much more frequent than either electron-electron 
or electron-ion collisions. Hence, in this case f(v) is the distribution function of the 
neutral background gas. 

The Lorentz gas model assumes elastic collisions with neutrals that are infinitely 
heavy and stationary, that is, 

F(v) = Ns(v) (4) 

where N is the number density of the neutrals. Equation (3) then becomes 


2 ) - 
^t/coll 



= N|v| y cr(x,|v |) f(v') 


-f(v) 


dST 


(5) 


where x is the angle between v and v ' and dfi' is the differential solid angle 
in v ' space. This collision model was used by Derfler (ref. 14) to derive the disper- 
sion relation for electrostatic waves in an unmagnetized plasma. Drummond, et al. 

(ref. 13), Suzuki (ref. 15), Idehara and Sugaya (ref. 16), and others have assumed that the 
collision cross section is isotropic, that is, ct is a function of |v| only. In this case 
equation (5) trivially becomes 



where v(v) = 47rNva(v) is a velocity-dependent collision frequency. 

Wilensky (ref. 18), in his work on electrostatic waves in an unmagnetized argon 
plasma, showed that the most important physics is in the energy dependence of the cross 
section and not in the angular dependence. He allowed for an angular -dependent cross 
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section and found that the wave dispersion differed very little from that given by the iso- 
tropic collision model. Consequently, ignoring the angular dependence as in equation (6) 
should be a very good approximation. In this report the collision models given by equa- 
tions (5) and (6) are used exclusively. 

In an attempt to put the Boltzmann equation in a more tractable form than the col- 
lision integral allows, a large number of approximate collision models have been used. 
Dougherty (ref. 19) in his investigation of the upper ionosphere used a model Fokker- 
Planck equation 

«„JL. 

\ 9t /coll 3v 

where v is a constant electron collision frequency and <f and T are a velocity and a 
temperature which Dougherty chose so that the generation-dissipation rates take on the 
proper values as discussed later. This collision model is more appropriately used when 
collisions result in only small changes in the velocity vector of the colliding particle and 
hence is most appropriately applied to coulomb collisions. This equation is derived in a 
more general form from equation (3) by Holt and Haskell (ref. 20) . On the other hand, 
Bhatnager, Gross, and Krook (ref. 21) replaced Boltzmann’s collision integral by a simple 
dissipation term 



which describes the relaxation of electron velocities to a local Maxwellian distribution 


(7-iT)f + i£4 

m dV 


(7) 


f 


m 



(9) 


Similarly, Allis, Buchsbaum, and Bers (ref. 22) assumed relaxation to some as yet 
unspecified equilibrium distribution function fg and wrote 



( 10 ) 


Note that the last three models all contain an average collision frequency v which must 
be specified, and equations (7) and (8) have other unspecified parameters. These unknown 
parameters should in practice be specified so that the principle of detailed balancing is 
obeyed. This procedure involves taking the following generation-dissipation rates of 
particles, momentum, and energy: 
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( 11 ) 



and choosing the unknown parameters so that these rates take on their proper values. 

For example, I should be chosen equal to zero if no particles are generated by the 
collision process. If the collisions are between particles of the same species, then the 
choice G = H = 0 should be made as Dougherty did, whereas for the Lorentz gas model 
described earlier H should be equal to zero but G not equal to zero. Achieving the 
correct generation-dissipation rates is not a problem with the collision models expressed 
by equations (3), (5), and (6) as they automatically take on the correct rates to all orders, 
and not just to the first three expressed by equations (11) to (13). 

Outline of the Present Analysis 

As mentioned earlier, the collision models expressed by equations (5) and (6) are 
used exclusively in this report. Waves proportional to exp (icot - ik • r ) are assumed 
and the Boltzmann equation with equation (6) is solved for the general wavenumber k . 

The solution is then specialized to k perpendicular to the magnetic field B, and the 
conductivity tensor and dispersion relations are obtained. The stability and dispersion 
of the longitudinal or electrostatic waves are investigated by numerically solving the 
appropriate dispersion relation. The stability character of the waves is determined by 
mapping the real wavenumber axis into the complex frequency plane by means of the dis- 
persion relation and then applying Derfler's frequency cusp criterion (refs. 23 to 26). 

In a somewhat simplified form sufficient for use here, if the mapping in the complex fre- 
quency plane is above the real frequency axis, the plasma is stable. If it goes below the 
axis, the plasma is at least convectively unstable. It is absolutely unstable only if a 
branch point can be located either on the mapping but below the real frequency axis or in 
a closed region above the mapping but below the real frequency axis. The Boltzmann 
equation with equation (5) is solved and the dispersion relations and conductivity tensor 
for k parallel to the magnetic field are obtained. The dispersion relation representing 
the transverse electromagnetic waves is numerically solved, and solutions of the disper- 
sion relations for various isotropic electron distribution functions are given. For propa- 
gation both parallel and perpendicular to the magnetic field, a monoenergetic electron 
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distribution is treated first and then a peaked energy distribution function having a peak 
of nonzero half -width is considered. Consideration is also given to a Maxwellian dis- 
tribution function and to a Maxwellian with a bump on the tail. It is assumed that elec- 
trons collide with nitrogen neutrals, but similar results would be obtained for any gas 
having a steep slope in its electron-neutral collision cross-section curve. In particular, 
the Ramsauer gases (argon, xenon, and krypton) would be appropriate because of their 
rapidly increasing collision frequency near the Ramsauer minimum. 

SYMBOLS 


A defined by equation (H5) 

A(v) defined by equations (27) 

Aj(^) defined by equation (37) 

a^ parameter used in equation (El 3) 

B defined by equation (H6) 

B magnetic field vector 

Bq,Bq zero-order magnetic field vector and magnitude 
Bj first-order magnetic field vector 

B^(x) Bernoulli polynomials (eq. (H9)) 

C Q number of electrons gaining energy 

Cj^ number of electrons losing energy 

C m constant, (eq. (78)) 

C s constant in peaked electron distribution (eqs. (73) and (78)) 

C(t) dummy function (eq. (B8)) 

c velocity of light 
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velocity parameter in equation (7) 


D,D 1 ,D 2 ,D e ,D ± dispersion relations 


D ’ D ij 

dfi' 


d 3 r 


defined by equation (29) 
element of solid angle in v space 
element of solid angle in v* space 
differential volume 


E 


E, 


E X ,Ey 


E, 


ei,e 


1’ 3 


f(v) 


electric field vector 

first-order electric field vector 

components of E^ along x and y directions 

symbolic for E x ± iEy and represents right- and left-hand polarized 
waves (eq. (122)) 

unsigned charge of electron 

unit vectors along X axis and Z axis 


e x ,ey,e z unit vectors along x, y, and z directions 


particle distribution function (eq. (3)) 


pFq^Fg generalized hypergeometric function 

F(z) dummy function (eq. (II)) 

F (<p) dummy function (eq. (E8)) 

f electron distribution function 


f H homogeneous part of electron distribution function (eq. (B5)) 
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f 

i, expansion variable of first-order electron distribution function (eq. (84)) 

f m Maxwellian electron distribution (eq. (9)) 

f g peaked electron distribution function 

fg zero- order electron distribution function 

f j first-order electron distribution function 

f components of first-order electron distribution function (eq. (M5)) 

f j transform of first-order distribution function 

G momentum generation rate (eq. (12)) 

G(v) defined by equations (36) and (52) 

g s defined by equation (97) 

g ± defined by equation (109) 

H energy generation rate (eq. (13)) 

Hj defined by equation (H10) 

h collisional frequency slope parameter used in equation (67) 

h(v) defined by equation (21) 

defined by equation (H7) 

I particle generation-dissipation rate (eq. (11)) 

I unit tensor 

I^(z) modified Bessel function 

*1’*2’*3 convenient symbols denoting specific integral used in appendix E 
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Im( ) imaginary part of variable 

i A 

j current density (eq. (33)) 

^ext driving current (eq. (121)) 

j x ,jy components of current density along X axis and Y axis (eq. (121)) 

jg( )d 2 ( )>• • *J n ( ) spherical Bessel functions of first kind 

j ± symbolic for j x ± ij y (eq. (122)) 

k ,k wavenumber vector and magnitude 

L,Ly defined by equation (30) 

l , m indices of summation 

m mass of electron 

N neutral particle density; electron number 

Nq number of electrons gaining energy 

Nj^ number of electrons losing energy 

n electron density 

n unit vector along B 

n i ,n j’ n k components of vector n 

P polarization vector 

P c probability of collision 

P^( ) Legendre polynomials 
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P z m () 


Re() 


v 0 

V 1 

v 2 

v 3 


associated Legendre polynomials 

pressure in equation (71); dummy variable in equation (H2) 

defined by equation (86) 

transfer collision cross section (eq. (87)) 

real part of variable 

space position vector 

symbolic notation given in equation (41) 

parameter in peaked electron distribution (eq. (73)) 

temperature 

defined by equation (91) 

time variable; dummy variable (eq. (B17)) 

dummy variable (eq. (B17)) 

velocity vector (eq. (3)) 

defined by equation (H9) 

velocity vectors 

component of velocity vector 

thermal speed; velocity at maximum of peaked distribution function 
thermal speed (eq. (78)) 

location of peak in peaked distribution (eq. (78)) 
component of velocity along Z axis 
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w 7 

x, Y,z 


x m-y m 


Y f,m<^> 

Z() 

z F () 

z i 

Of 

r() 

r 

5(v) 

^r,s 

e ijk 

e 0 

V 

e,e' 


weights of Gauss quadrature rule (eq. (G2)) 

weights of Gauss quadrature rule along contour (eq. (18)) 

coordinate axes 

nodes of Gauss quadrature rule (eq. (G2)) 
parameters in complex contour (eqs. (12) and (13)) 
spherical harmonic 
defined by equation (N8) 

Fried function 

nodes of Gauss quadrature rule along contour (eq. (17)) 
proportion of electron distribution function that is Maxwellian (eq. (78)) 
jOfjQ collision frequency function parameters (eq. (70)) 

defined by equation (98) 
gamma function 

parameter in complex contour (eq. (13)) 

Dirac delta function 

Kronecker delta 

permutation tensor 

free- space dielectric constant 

normalized frequency, u> v j<^ c 

spherical coordinate angles in velocity space 
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k Boltzmann’s constant 

\ defined by equation (55) 

\ n ^ s defined by equation (99) 

Xq defined after equation (N4) 

\ defined by equation (117) 

jjLq permeability of free space 

v collision frequency (eq. (6)) 

transfer collision frequency (eq. (88)) 
collision frequency at infinite electron velocity 
(i/> average collision frequency (eq. (77)) 

| dummy variable (eq. (59)) 

77,77^ polarization tensor (eq. (Cll)) and its components 

p charge density 

a(x>v) collision cross section 

a n (v) expansion constant of collision cross section (eq. (J2)) 

o^CTy conductivity tensor 

dummy variables (appendix B) 

cp spherical coordinate angle in velocity space; also a dummy integration 

variable 

^(<P>v) defined by equation (28) 


13 



X 


scattering angle 


oo frequency angle 

3 c> w c 

W H 


Subscript: 

coll collision 

Vector symbols without arrows denote the magnitude. An asterisk denotes a com- 
plex conjugate and a circumflex denotes a transform. 

THEORY OF WAVES PROPAGATING PERPENDICULAR 
TO THE MAGNETIC FIELD 

In this section dispersion relations for waves propagating perpendicular to an 
applied magnetic field Bq are derived. A Lorentz gas model is assumed, and the 
electron-neutral collision cross section is assumed to be isotropic. The first order dis- 
tribution function is derived and from it the conductivity tensor for the general case of 
IT at any angle to Bq is obtained. This tensor is then specialized to propagation per- 
pendicular to the magnetic field, and the dispersion relations are obtained. A discussion 
of some of the analytic properties of the dispersion relation that describes electrostatic 
or longitudinal waves is given. 

Conductivity Tensor for General k With Isotropic Collisions 

In the Introduction it was stated that the Boltzmann equation for electrons colliding 
elastically and isotropically with infinitely heavy neutrals is 

— + v • -4S ” — (e + v X b) • -^4 = -^(v)f(v) + C f(v')dS2' (14) 

8t dr ™ ' av 4 tt J 


cyclotron frequency (eq. (22)) 
hybrid frequency, \J^ + w p 


1/2 


plasma frequency, ^ne 2 yW 0 ^ 
abbreviated notation for oo - iv (eq. (23)) 
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where dS2' is an element of solid angle in v ' space. This equation is linearized by 
making the "ansatz" (initial assumption) 


f(v) = f Q (v) + fj(v) 

E = 0 + Ej \ 


(15) 


and in the usual manner consider the first-order terms (subscript 1) to be small pertur- 
bations about the much larger zero-order terms (subscript 0). Then obtain the zero- 
order equation 


- ^ x ®o) ■ $ - -•*<> + h r v™* 


(16) 


and the first-order equation 

afj „ , _ x di 

dr av 


af 

"at 


3f 3f c* 

- + v . — - - — (v X B„\ • — - = — (e. + v X • — - - id:, + — \ f, (v ')dfT 

: mV 0) mV 1 11 dv 1 4tt J 1 


(17) 


In appendix A the solution of the zero-order equation is found to be any isotropic distri- 
bution function fg(v) and hence the term (v x BjJ • di^jd^ in the first-order equa- 
tion is identically zero. To solve equation (17), use is made of a Fourier transform in 
space defined by 


f (v ,TE ,t) = y f(v",?,t)e i ^ * r d^r 
and a Laplace transform in time defined by 


(18) 



(19) 


where s is chosen equal to iw. Applying these transforms to equation (17) or equiva- 
lently assuming that all the variables are proportional to exp(ici>t - ik • r ) yields 


af. 


i(w„ - k • v)fj - (v x u> c ) • — - = h(v) 


(20) 
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where 



( 21 ) 


(22) 


00 ^ - Cl) - il'(v) 


(23) 


The Laplace transform analysis was used so that the rules of the inverse Laplace 
transform could be used later to define the dispersion relations properly. Even though 
no initial conditions are allowed for, the dispersion relations that are derived will 
properly describe any waves that are excited and will predict the onset of wave 
instabilities. 

Equation (20) is solved by the method of characteristics in a manner similar to 
that used by Oraura (ref. 27). A parameter r is introduced so that v = v"(t). The 
characteristic equations then become 


d|ti=-(vx 3c ) 


/ - df 1 
i( w v - k • vjfj + — = h(v) 


(24) 


(25) 


In appendix B these equations are solved and the solution is found to be 


f 1 = — 
1 mov 


~ — f A(v ' )dfi' 

, 47TCO r J Jn 

A(V) + ^ 2 


47TCU. 


jy dS2 , 


(26) 


where 
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( 27 ) 
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( 28 ) 


= i <p(v -jy) -E • LM : v 

W C 

Hi] = D ij = n i n i + ( 6 ij - n i n j) cos p + e ijk n k sin * 

Hij = L ij = n i n ]^ + ( 6 ij ■ Vi) sin ^ - e ijk n k (cos * - 1} 

_ e - B 

= — B n = a) — = gj n 
c m 0 c Bq c 


(29) 

(30) 


(31) 


By direct substitution of equation (26) into equation (20), one can show that the solution 
has been obtained. 

Equation (26) was also obtained by Idehara and Sugaya (ref. 16) in a somewhat more 
complicated form. However, they were unable to simplify it to obtain relatively simple 
dispersion relations of wide validity, as will be done here. 

In the collisionless limit equation (26) yields 


U = -4- A(v) (32) 

1 mu> c 

This result was obtained by Omura (ref. 27) and others and can be used to derive 
Bernstein's (refs. 28 and 29) dispersion relation describing the propagation of cyclotron 
waves in the Vlasov limit. 

In appendix C, Maxwell's equations are stated and the general dispersion relation 
is derived for a conductivity tensor a defined by 


j = a • E = -en ^ v f ^(v)d^v 

Using equation (33) and taking B = Beg in equation (26) yields 


(33) 


a ir 


Cl) 

“a> 


2 

5 T fo Vf ‘ )<V> F V i A j <ir H d ' r 


"p e 


4770) 


/-fc°0 

5J 0 vl ' Wf o< v > 


^ Aj (v)dfi y y 0 clO 


1 + OM G(v) 


dv 


(34) 
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where 


"n = — 

P me„ 

G(v) . f f“ e-^’^dv d n 

47rw c 

v-> =1“ e-* ( ^v k » 

(poo - k • L • v 
H<P,V) = i ^—r: 


L = 


sin cp 

-(1 - cos cp) 
0 


(1 - cos cp) 0 

sin (p 0 

0 <p 


D = 


COS <p 

-sin cp 


sin (p 
cos (p 


0 0 


0 

0 

1 


(35) 

(36) 

(37) 

(38) 

(39) 


This tensor is the exact conductivity tensor appropriate to the model assumed. In the 
next section this tensor is put in a usable form after assuming k _L B, and the wave 
dispersion relations are then obtained from Maxwell's equations. In a later section 
k || B is assumed and again the simplified conductivity tensor and wave dispersion 
relations are obtained. 

Conductivity Tensor and Dispersion Relations for k J_ B 

In this section the conductivity tensor given by equation (34) is simplified by con- 
sidering propagation strictly perpendicular to the magnetic field, for example, k = ke^. 
For convenience, the symbolic notation 
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(41) 


oo 

W " I 


( kv /"c) 


2 Z 


1^1 Of 2 - l) . . . U 2 - * 2 ) 


where 77 = <i> v j<o c will be used. 

In appendixes D and E after much algebraic manipulation, the 
is given by 


where 
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conductivity tensor 


(42) 


(43) 


(44) 


(45) 


(46) 
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Equations (42) to (46) give the conductivity tensor in a form suitable for computations. 
Shimomura and Mitani (ref. 17) attempted to derive this tensor also, but they made an 
error in solving the first-order Boltzmann equation given here by equation (20). Their 
error was of such a nature, however, that they still obtained the correct expression for 
Ogg. It agrees with equation (46) when allowances are made for differences in notation. 

To obtain the dispersion relations, assume k = ke^ in equation (C17), and use 
equation (42) to obtain 



D 2 (o>,k) = ll + 1 + — - l—f + (— 


The waves described by D^(w,k) are characterized by || Bq and are usually called 
ordinary waves, since they are unaffected by the magnetic field in the cold plasma limit. 
This same dispersion relation was derived in an approximate form and discussed by 
Idehara and Sugaya (ref. 16) for a monoenergetic electron distribution. Shimomura and 
Mitani (ref. 17) also obtained equations (48) and (46) and found that a monoenergetic dis- 
tribution function leads to collisional instabilities if the electron-neutral collision fre- 
quency increases rapidly enough with velocity. This dispersion relation is not discussed 
further in this report. 

In the limit of very small phase velocity, that is, kc/ co »> 1, equation (49) yields 

D 2 MM + ^- = 0 (50) 

0 
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Using equation (43) and dropping the subscript 2 on D then yields 


where 


D(co,k) = 1 + 47T 



v f o"„ G W 

c y + ii^G(v) 


dv = 0 


G(v) = £ 


(kv/u^) 


2 1 


=! (& + 1)(?7 2 - 1) ... (?7 2 - Z 2 ) 


( 51 ) 


(52) 


The waves described by this equation are characterized by | k and in the collision- 
less limit have been variously called extraordinary, longitudinal, electrostatic, cyclotron, 
or Bernstein waves. They were discussed in a review paper by Crawford (ref. 30) in 
1967. The dispersion relation describing these waves is generally derived in a much 
simpler manner by using the electrostatic approximation ab initio. This derivation has 
been made in appendix F where equation (51) is obtained more directly. This dispersion 
relation is discussed further in the next section and is later solved for several inter- 
esting cases. 


Discussion of Electrostatic Dispersion Relation 

In the preceding section, the dispersion relations for waves propagating perpendic- 
ular to an applied magnetic field in a Lorentz gas plasma were derived. These disper- 
sion relations were exact in the framework of the collision model assumed, that is, an 
isotropic velocity-dependent electron-neutral collision cross section. The quasi-static 
approximation was then made to obtain the dispersion relation for cyclotron or 
Bernstein waves. The properties of this dispersion relation will now be discussed. 

From the series expansion given by equation (52), G(v) can be written alterna- 
tively as a generalized hypergeometric function 


G(v) = -1 + 


2 F 3 


1 , 


1/2; 3/2, 1 - 77, 1 + 77; 



(53) 


Consequently, the theory of asymptotic expansions of hypergeometric functions can be 
used to help evaluate G(v) (ref. 31). 


In the limit of zero magnetic field, that is, w_ — 0, equation (51) reduces to 


/ w n\ 2 r oo vfj, (l 


D <“’ k > = 1 - 4 'i?) 


A coth'^A 


1 + — coth"^A 
kv 


dv = 0 


(54) 
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nil mi 


where 


X 



(55) 


This dispersion relation was obtained by Drummond et al. (ref. 13) and Derfler (ref. 14) 
and represents electrostatic waves in an isotropic Lorentz gas plasma. Hence, equa- 
tion (51) has the correct form in the limit of a vanishing magnetic field. 

On the other hand, in the presence of a magnetic field, but in the limit of zero col- 
lisions, equation (51) yields 

fu> \ 2 p°° 

D(o;,k) = 1 + j vfgG(v) dv = 0 (56) 


If fg is taken to be Maxwellian, that is, 



then equation (56) yields 


where 


D(w,k) 


V \ 2 y 1 -3-5 . . ■ (a - 1) 

kv o/ 2 l (r] 2 - l) . . . (?7 2 - l 2 ) 


(57) 


(58) 


V = 


00 - IV 


00 , 


This relation can be shown to be equivalent to Bernstein's dispersion relation (refs. 28 
and 29), given by 


where 




2 


(59) 
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However in practice, equation (58) has proven to be much more efficient for computation. 
The collisionless electrostatic dispersion relation, given by either equation (58) or 
equation (59), is characterized by undamped pass bands at frequencies just above each 
electron cyclotron harmonic separated by bands that are very heavily damped. In the 
limit of zero magnetic field, neither of these equations predict Landau damping since the 
infinite series in both equations represents only the contribution from an integral along 
the real axis in velocity space. However, in this limit the electrostatic waves are 
damped as Landau (ref. 32) and others (ref. 33) have shown. To recover Landau damping 
properly it is necessary to indent the velocity contour around a singularity that crosses 
the real axis as the magnetic field vanishes. Baldwin and Rowlands (ref. 34) have shown 
how Landau damping can be recovered by rewriting equation (59) in an integral form and 
eventually obtain the Landau dispersion relation 


D(u>,k) = 1 - 



(60) 


where Z' is the derivative of the Hilbert transform of the Gaussian as tabulated by 
Fried and Conte (ref. 35). However, the theory given here quite easily yields equa- 
tion (60) from either equation (56) as the magnetic field vanishes or from equation (54) 
as collisions vanish, if it is assumed, of course, that fg is chosen to be Maxwellian. 
When doing this, one must keep in mind that the dispersion relations are defined for k 
real, and co on a Laplace integral path. For other k and u> values, the appropriate 
analytic continuation must be taken by choosing the correct contour in velocity space. 

The dispersion relation given by equation (51) is considerably different from those 
given by cruder collision models. Use of the simple constant collision frequency model 
of Allis et al. (ref. 22), that is, 


(Dcoir^-*) < ei > 

has the effect of replacing u> by u> - iu in the collisionless equation. In an attempt to 
improve on this collision model, a variation of the Bhatnagar-Gross-Krook (BGK) colli- 
sion model has been used (ref. 36), 



(62) 


where ng is the zero-order electron density and n^ is the first-order electron den- 
sity. In addition to the replacement of w by oo - iu, this model also replaces 
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by [l - i (I'/ctOjoip in the collisionless equation. Note that neither of these collision 
models correctly produces the denominator in the integrand of equation (51). This 
denominator is of major importance in determining the characteristics of the solutions 
of the dispersion relation as it tends to cancel the resonance effect of the integrand when 
the frequency is near a multiple of the cyclotron frequency. Hence one must conclude 
that the collision models expressed by equations (61) and (62) can only be applied in the 
near-collisionless and off-resonance limit. 

Next note that equation (51) contains only even powers of k, and hence the following 
symmetry relation holds: 

D(w,k) = D(o>,-k) (63) 

Consequently, the roots of D(co,0) = 0 are branch points in the mapping of k into the 
complex frequency plane. The location of these points is of prime importance in deter- 
mining the stability character of the waves (refs. 23 to 26). From equation (51) also 
obtain the symmetry relation 


D(w,k) = D*(-ci>*,-k*) (64) 

This relation is, of course, a statement of the Kramers-Kronig (ref. 37) relationship and 
is a direct consequence of the causality principle which was built into this theory by 
using the Laplace transform analysis in time. 

Combining equations (63) and (64) yields then for k real 

D(w,k) = D*(-w*,k) (65) 

Hence, only complex frequency solutions need to be determined in, for example, the right 
half -plane. These properties are important when the dispersion relation (eq. (51)) is 
solved for several different electron distribution functions in the next section. 

SOLUTIONS OF THE ELECTROSTATIC DISPERSION RELATION FOR 
PROPAGATION PERPENDICULAR TO THE MAGNETIC FIELD 

In this section the dispersion relation is solved for electrostatic waves propagating 
perpendicular to the magnetic field for several isotropic electron distribution functions. 

A series of new collisional modes are found that depend for their existence on both a 
highly non-Maxwe Ilian electron distribution function and a velocity-dependent collision 
frequency. Both convective and absolute instabilities are discovered when the electron 
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distribution function is monoenergetic. A Lorentz gas plasma with nitrogen neutrals is 
also considered and instabilities are again found when a peaked electron distribution 
function is assumed. When a Maxwellian electron distribution function is chosen, it is 
found that the wave dispersion predicted by the isotropic collision frequency model is 
very close to that predicted by a constant- collision-frequency BGK model. Also an 
electron distribution function consisting of both peaked and Maxwellian components is 
briefly considered and the stabilizing influence of the Maxwellian population is 
demonstrated. 


Monoenergetic Electron Distribution 

For the initial analysis of the dispersion relation (eq. (51)) it is assumed that the 
zero-order electron distribution function is monoenergetic, that is, 


f n( v ) 


_ 6 ( v - v o) 




(66) 


Substituting this relation into equation (51) leads to the dispersion relation 


’ a;? \ ^ fa) ^(w - if)^G + VqG') + iuG(G - hG - hjj + v^G^I 
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The series representations given for G and G’ were found to converge very rapidly; 
thus, equation (67) was easily solvable by numerical methods. 

Idehara and Sugaya (ref. 16) considered the same problem but they did not derive 
the relatively simple dispersion relation given by equation (67). They instead made 
approximations that are valid only for very small wavenumbers (or large wavelengths). 

Before solving equation (67) the roots of D(ct>,0) = 0 are investigated, since from 
a preceding section, these roots are known to correspond to branch points in the complex 
frequency plane. By expanding G and G' in equation (67) and keeping only terms to 
lowest order in k, the dispersion relation yields 


+ ufi (-4^) + w*^-6i^ - 2co 


+ w 


(z , 2 + a>2) 2 + ( 3 „ 2 


+ w; 



= 0 


( 68 ) 


This equation, being a complex fifth-order polynomial in u>, has five complex roots. 
However, because of the symmetry relations expressed by equation (65), there can be 
only three independent solutions. One root must have a zero real part, and each pair of 
the other four must be related by = -u>g where and denote one pair of 

the four remaining roots. One of these pairs has a real part near the hybrid frequency 

o> H = \J^ + a>p and the other has a real part near the cyclotron frequency. The precise 

location of these roots will, of course, depend on the plasma parameters chosen. If 
these roots lie below the real axis, then an absolute instability is indicated, whereas if 
they are above the axis, there is no absolute instability but there may or may not be a 
convective instability (refs. 23 to 26). By solving equation (68) for various values of the 
slope parameter h, it is found that the root with the real part near the cyclotron fre- 
quency leads to an absolute instability for large positive values of h, for example, 
h > 3. However, for large negative values of h, the roots with real parts zero and near 
the hybrid frequency both lead to an absolute instability. In addition to these solutions, 
one can show analytically that as k — 0, equation (67) is solved by w = 0 and 
w = ±nw c + iv where n = 2,3,4,. . . . 

With the help of these solutions, it is possible to numerically solve the dispersion 
relation (eq. (67)) with k real and greater than zero to obtain the desired complex fre- 
quency solutions in the right half-plane. This procedure leads to the mappings of the 
real wavenumber axis into the complex frequency plane shown in figures 1 to 4. 

In figure 1 for h = 4, an absolute instability at Re(co) ~ o> c is predicted since one 
of the branch points (indicated by the small circles) is in the lower half -plane. In this 
figure only the first three terms of an infinite series of modes are plotted. The modes 
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Figure 1.- Real wavenumber axis mapped into complex frequency plane 
with a positive slope parameter. The branch point (indicated by 
a small circle) below the real od axis indicates an absolute 
instability. The upper set of modes are the modified collisionless 
modes. The bottom set are the colli si onal modes. The hybrid fre- 
quency is indicated by the solid square symbol on the real <x> axis. 

near higher cyclotron harmonics are similar to the two near oo = 3c o Q in figure 1. The 
two modes with Re(cu) = 0 are collision-dependent modes and in the collisionless limit 
they both coalesce to the point cu = 0. They are not of much interest here since in a real 
plasma they would be masked by ion effects. The three upper modes in figure 1 are 
modified collisionless modes. In the collisionless limit, these modes have co real solu- 
tions and would hence map into lines on the real frequency axis. The three lower modes 
in this figure are collision-dependent modes, and in the collisionless limit they coalesce 
into the points c o = ±nco c (n = 1,2,3,. . .). Note particularly the absolutely unstable mode 
with Re(cu) ~ cu c . The location of the branch point associated with this mode is given by 
one of the roots of the fifth-order polynomial (eq. (68)). If either v = 0 or h = 0, the 
order of this equation is reduced by two as can be seen by reobtaining the new equation 
from equation (67). The roots that vanish have real parts near ±co c and hence it is 
concluded that these modes are introduced by collisions. The branch points associated 
with the other collisional modes (and with the set of modified collisionless modes also) 
are located at co = ±nco c + iv (n = 2,3,4,. . .). One can show analytically that if either 
h or v goes to zero, the collisional modes vanish and the branch points are then 
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Figure 2.- Real wavenumber axis mapped into complex frequency plane 
with a negative slope parameter. The branch point at the hybrid 
frequency is in the lower complex frequency plane so that the 
plasma is absolutely unstable. Rote that the modified collision- 
less modes are below the colli si onal modes rather than above as 
in figure 1. Hybrid frequency is indicated by a solid square 
symbol. Small circles denote branch points. 

associated only with the collisionless modes. These branch points cannot lead to abso- 
lute instabilities for any set of plasma parameters. At most, these modes may be con- 
vectively unstable unless there is another branch point below the real frequency axis but 
above the mapping of the real k. axis. No branch points of this type were found. Any 
other branch points in the lower frequency plane may be shown to be ignorable. For 
large k, the mappings appear to approach the points c o = nc o c + \v (n = 1,2,3,. . .), but 
this result was not shown conclusively. It is not expected, however, that anything of 
physical interest occurs at k values larger than those shown. 

In figure 2 it is seen that for h = -4, there are two branch points in the lower half 
plane. Note that the position of the modes has reversed. The modified collisionless 
modes are on the bottom, and the collisional modes are on the top. The unstable modes 
are hence the modified collisionless modes, and there is an absolute instability at the 
hybrid frequency. The modes near Re(co) = 0 are also shown here, but as before, are 
of only minor interest. Henceforth, they will not be considered. 
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Figure 3«- Real wavenumber axis mapped into the complex frequency plane. 

The branch points are all above the real to axis but the mapping 
does dip below this axis and indicates a convective instability. 

Hybrid frequency denoted by a solid square symbol; branch points are 
denoted by small circles. 

In figures 3 and 4, cases that are convectively unstable are shown. Note that the 
branch points all lie above the real frequency axis so that there is no absolute instability. 
However, the mapping of the real wavenumber axis dips below the real frequency axis and 
indicates a convective instability. The unstable modes in both of these figures are colli- 
sion dependent and vanish in the collisionless limit where the waves supported by the 
plasma are known to be stable (ref. 36). Furthermore, for large collision frequencies a 
damping effect is evident as shown by the fact that the v = 0.05w c case in figure 4 is 
more unstable than the v = 0.1w c case in figure 3. Hence there must be some inter- 
mediate collision frequency (and hence neutral gas density) for which the growth rate 
Im(o>) would reach a maximum. Note also that figures 3 and 4 are for h = 2. Values 
of the slope parameter in this range are common in gases that are used in laboratory 
discharge plasmas. For this reason it appears that this theory will find important 
applications in predicting cyclotron wave propagation through plasmas. More realistic 
electron distributions are treated in subsequent sections. 

It should be pointed out at this point that a spherical shell distribution as considered 
here is unstable even in the collisionless limit for some ranges of the plasma parameters 
(ref. 36). However, for the cases investigated here and in the remaining part of this sec- 
tion, the collisionless limit predicts only stable or evanescent waves. Hence, one can 
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Figure 4.- Real wavenumber axis mapped into the complex frequency plane 
for a convectively unstable case. Hybrid frequency indicated by a 
solid square symbol; branch points are denoted by small circles. 


truly say that the electron-neutral collision process is responsible for the predicted 
instabilities. 


Peaked Electron Distribution 

In a previous section equation (51) is solved for a monoenergetic electron distribu- 
tion function and collision-induced instabilities are found for some ranges of the slope 
parameter 


v d v_ 
v dv 


(69) 


In particular, unstable modes near the cyclotron harmonics are found for values of h 
near 2 and above. For h negative, unstable modes with real frequency near zero and 
near the hybrid frequency are found. Hence it seems likely that collisional instabilities 
of this type should be found in a plasma if there were present (1) a background of neutrals 
for which h > 2 over some region or regions of velocity space, and (2) a sharply peaked 
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electron distribution function centered in one of these regions. In the remainder of this 
section it is shown that this is indeed possible. 

The measured electron-neutral collision frequency for nitrogen exhibits a sharp 
peak (refs. 38 and 39) as shown in figure 5. The parameter h ranges up to 5; thus 
criterion (1) is met. To use the experimental values in the dispersion relation, first fit 
a rational polynomial 



0 12 3 4 


VELOCITY IN JesT 

Figure 5 *- Rational function fit to experimental electron-neutral 
collision frequency for nitrogen (ref. 38) and slope parameter 
from rational function. 
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to the experimental probability of collision for nitrogen measured by Brode (ref. 38). 
Then obtain the collision frequency (ref. 39) from 

„ W ,™pP c v (71) 


where p and T are, respectively, the pressure in torrs and the temperature in K of 
the nitrogen. Note that since p is proportional to T, only the neutral particle number 
density is a parameter in the collision frequency function. Instead of number density, 
however, the collision frequency at infinite electron velocity defined by 

y vP 

v(v) = - m - c (72) 


a 


1 


will be used. A least-squares fit of equation (70) to the experimental probability of col- 
lision for nitrogen (and mercury) (refs. 38 and 40) led to the parameters given in table I. 
The quality of the fit for nitrogen can be seen in figure 5. 


TABLE I.- LEAST SQUARES FIT TO EXPERIMENTAL PROBABILITY 
OF COLLISION FOR NITROGEN AND MERCURY* 

= a (v 2 + o- 2 v + <* 4 ) (v 2 + « 6 v + Qfe) 

1 (v 2 + Q(3V + a 5 ^v2 + a7 v + ^(v-Q !l0 ) 
where P c is in cm 2 and v is in \/eV 



Nitrogen 

Mercury 

a l 

274.304 

381.030 

“2 

-2.58343 

.531224 

a 3 

-3.16090 

-.107263 

“4 

4.28551 

.516368 

“5 

2.55450 

.475607 

a 6 

-3.35392 

-.717384 

a 7 

-5.25412 

-.840921 


3.09938 

.132412 


12.5096 

.180940 

“lO 

-4.55542 

-7.01816 


^References 38 to 40. 
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To meet criterion (2), the distribution function, 


f 0 (v) = f s (v) 



( 73 ) 


where 


^(s + 3)/2 



is used. This distribution function has a peak at Vq and is shown plotted in figure 6. 
The half-width of the peak is a transcendental function of s and is not conveniently used 
as a parameter. However, the half -width increases (1) proportionally with Vq and/or 
(2) as s becomes smaller. The electron distribution function (eq. (73)) is normalized 
according to 

y f 0 (v)d 3 v = 1 (74) 


Once the distribution function fg(v) and the collision frequency function u(v) 
are specified, the dispersion relation in equation (51) can be solved. It is necessary to 
evaluate the velocity integral numerically but this is easily done by using a Gauss quadra- 
ture rule. The details of this integration procedure are given in appendix G. When 



Figure 6.- Normalized, peaked distribution function. 
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carrying out this integration, it was necessary at times to use an asymptotic expansion 
for G(v). Since G(v) was known to be a hypergeometric function (see eq. (53)), it was 
possible to utilize the asymptotic expansions of generalized hypergeometric functions 
given by Meijer (ref. 31). This expansion is given in appendix H. Note also that the 
integrand in equation (51) will have poles wherever the denominator is zero. Once the 
poles are located in velocity space, it is necessary to determine the correct integration 
contour before trying to solve the dispersion relation. This determination is made by 
recalling that equation (51) is defined for k real and co on a Laplace integral path 
(LIP). Usually, the real velocity axis proved to be a valid contour, but occasionally, it 
was necessary to integrate along a deformed contour as described in appendix I. The 
scheme described worked very well but is only one of many that could be used. 

To determine the stability of the waves predicted by equation (51) for the case being 
considered in this section, it was decided to map the real wavenumber axis into the 



(a) = 0.5a> c . 

Figure Mapping of real wavenumber k into the complex frequency 
plane. Plot shows absolutely unstable waves. Hybrid frequency is 
indicated by a solid square symbol; branch points are denoted by 
small circles. 
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Re (u/u c ) 

(t>) v ro = 0.1a> c . 

Figure 7*- Concluded. 

complex frequency plane. If instabilities of an unknown nature are present, this method 
seems to be a more useful approach than, for example, mapping from the frequency plane 
into the wavenumber plane. 

Figures 7 to 9 show the results of such a mapping. Because of the symmetry 
relations expressed by equations (63) and (65), it is necessary only to consider k > 0 
and to look for roots in the right half frequency plane. Figure 7(a) shows a case that is 
absolutely unstable because of the branch point at co « cu c . The electron distribution 
function peak is at Vq = 1.25\/eV and from figure 5 it can be seen that this point is the 
approximate location of the most positive slope of the collision frequency function. The 
three lower modes in figure 7(a) represent the new collisional modes described earlier 
and vanish as the collisionless limit is approached. The three upper modes are modifi- 
cations of the collisionless modes but are damped here unlike in the collisionless limit. 

It is interesting to note that the modes shown in figure 7(a) all come from their own 
branch points at k = 0, whereas in the monoenergetic electron distribution case, there 
were sometimes two modes arising from a single branch point. 

The case shown in figure 7(b) is identical to that shown in figure 7(a) except for a 
lower collision frequency. Note that there are now more unstable modes; thus, it 
appears that high collision frequencies tend to damp even the collisional modes. Since 
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Figure 8.- Mapping of real wavenumber k into the complex frequency 
plane. Plot shows convectively unstable waves. Hybrid frequency 
is indicated by a solid square symbol; branch points are denoted 
by small circles. 

these modes also vanish at zero collision frequency, there is a particular collision fre- 
quency, and hence a particular neutral gas pressure, at which the collisional cyclotron 
waves are most unstable. 

In figure 8 is shown a case where the electron distribution function half-width is 
approximately twice that shown in figure 7(b). (See fig. 6.) Note that the first mode has 
changed from absolutely unstable to convectively unstable and the third mode has become 
stable. Hence, electron distribution functions with wider peaks tend to be more stable 
than those with narrow peaks. 

In figure 9 a case with the electron distribution function peak located on the negative 
slope of the collision frequency curve is considered. As in the case of the monoenergetic 
distribution function, the location of the modes is reversed and the collisional modes are 
very heavily damped. However, the case in figure 9 is a stable situation unlike the one 
presented in figure 2. Consequently, if the peaks were made sufficiently narrow, the case 
given in figure 9 would become absolutely unstable because of the branch point with real 
frequency near the hybrid frequency. However, the case of a positive slope parameter 
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Figure 9 >- Mapping of real wavenumber k into the complex frequency 
plane. Plot shows stable waves. Hybrid frequency is indicated 
by a solid square symbol; branch points are denoted by small circles. 


seems to be of more physical interest because of the large number of gases having 
regions of positive slope in their collision frequency curves. Consequently, the negative 
slope case will not be investigated further here. A case with large half-width, s = 1 
and Vq = 1.05\/eV, was also considered briefly. It was found that the collisional mode 
with real frequency near the cyclotron frequency was convectively unstable but that all 
other modes were stable. 

Earlier it has been shown that collisional instabilities are possible in a plasma if 
a large percentage of the electrons are in a region in velocity space where the collision 
frequency increases (or decreases) rapidly with velocity. If the collision frequency 
slope is positive, there may be an absolute instability at the cyclotron frequency and con- 
vective instabilities at the higher harmonics. If the collision frequency slope is negative, 
there may be an absolute instability at the hybrid frequency and convective instabilities 
at frequencies between the cyclotron harmonics. The occurrence and strength of these 
instabilities depend very strongly on (1) the width of the peak in the electron distribution 
function and (2) the slope of the collision frequency in the region of the peak. It is diffi- 
cult to give more precise statements about the occurrence of these collisional instabili- 
ties; thus it is necessary to consider each case separately. 
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Maxwellian Electron Distribution 


In this section equation (51) is solved with the same collision frequency function as 
in the last section, that is, equation (72), but instead of a peaked electron energy distri- 
bution function a Maxwellian distribution function 


f 0 ( y ) 



(75) 


is used. 

In figure 10 the real wavenumber axis is mapped into the complex frequency plane 
and the modes shown are of the modified collisionless type. These modes cannot be 
unstable since a Maxwellian plasma must be stable. As k becomes large, the mappings 
apparently approach the points u> = nu> c (n = 1,2,3,. . .). This statement, however, was 
not shown conclusively either by numerical or analytical means. The collisional modes, 
that were often unstable for the peaked electron distribution function, were not found in a 
Maxwellian plasma. Apparently, these collisional modes are present only if the electron 



Figure 10.- Mapping of real wavenumber k into the complex frequency 
plane for a Maxwellian electron distribution. Hybrid frequency is 
denoted by a solid square symbol; branch points are denoted by 
small circles. 
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distribution function has a region of positive slope. To investigate the disappearance of 
these modes further, the isotropic distribution function given by 


V v > = 



(v < v Q ) 
(v > v 0 ) 


(76) 


was briefly considered and it was found analytically that the collisional mode did not 
exist. In support of this argument, Bekefi et al. (ref. 41) found that, in general, there 
must be a region where afg/av > 0 for collisional instabilities to exist. However, they 
used a radiation temperature approach and could not say anything about the existence or 
nonexistence of the actual collisional mode. 

Since the waves present in a Maxwellian plasma are stable, the complex k roots 
for real co may be interpreted, according to established practice, as spatially damped 
or growing waves. This representation is relevant to laboratory experiments since one 
often excites the plasma at some real frequency and measures the complex and real parts 
of the wavenumber. The results of such a mapping are shown in figure 11. The case 
shown here was solved also by Tataronis (ref. 36) by using a constant collision frequency 
BGK model, that is, equation (62). Because of the presence of the velocity-dependent 
collision frequency in equation (51), the temperature of the Maxwellian electron distribu- 
tion. function did not normalize out as in his case, and so a thermal speed of Vq = 1.5\/eV 
was chosen as representative of values found in a laboratory plasma. For the collision 
frequency parameter v ^ a value was chosen so that the average collision frequency 
given by 

(v) - ^ v(v) f Q (v) 47rv 2 dv (77) 

was the same as the constant collision frequency used by Tataronis. From figure 11 it is 
seen that the two models differ but not substantially. There appear to be no important 
differences, unlike in the case of the peaked distribution function where the collisional 
modes do not exist for other collision models. Hence, it would seem that for a 
Maxwellian plasma, it is not necessary to use the exact collisional model considered 
here unless very precise results are desired. 


Sum of Maxwellian and Peaked Electron Distributions 

Since it was found that a peaked distribution function could lead to collision-induced 
instabilities, one might wonder what portion of the distribution function must be 
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Maxwellian before the waves become stable. To answer this question, the electron dis- 
tribution function is chosen to be 


f 0 (v) = “C m e 



+ (1 - a)C s 




(78) 


where 


C m = (^ v i) 


-3 


c s - 


_ (f) 


(s+3)/2 


2ttv 2 rf s + 3 


Note that a gives the proportion of the total electron population that is Maxwellian. 
This distribution function also obeys the normalization criterion (eq. (74)). 



Figure 12.- Stability plot for Maxwellian plus peaked distribution 
function. Mode with Re(co) ~ 




Figure 13.- Normalized distribution function consisting of peaked 
and Maxwellian components (eq. ( 78 )). Points A and B in 
figure 12. 


Parameters were then chosen that were known to be absolutely unstable, that is, 
s = 30, v 2 = l^Sy^eV, and ^ = 0.1o> c . The dispersion relation was solved for various 
values of a and v^ to obtain the instability boundaries in the a,v^ plane shown in 
figure 12. The mode considered here was the collisional mode that has Re(cu) ~ cn c . 

On the lower curve in figure 12, k = 0 since the branch point in the complex frequency 
plane is given by the root of D(cu,0) = 0. On the upper curve the mapping of the real k 
axis into the frequency plane dips down and just barely touches the real frequency axis 
but does not go below it. At the point of contact it was found that k ~ 2.5cu c /v2. 

Since an absolute instability occurs at zero group velocity (dco/dk = 0), it is not too 
surprising to find that low -temperature Maxwellians affect the waves a great deal since 
there would be proportionally more electrons near zero velocity. However, at higher 
temperatures it is somewhat surprising to see that as much as 95 percent Maxwellian is 
required to stabilize the unstable collisional mode. In figure 13 is shown a plot of the 
distribution function (eq. (78)) for a case on each of the boundary lines in figure 12 
(points A and B). Even though there are appreciable numbers of Maxwellian electrons, 
the distribution function is still highly non-Maxwellian. 
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Synopsis of Numerical Results 

Briefly, it has been found that a rigorous treatment of electron-neutral collisions 
in a non-Maxwellian Lorentz magnetoplasma predicts a new series of electrostatic modes 
that vanish both in the collisionless and Maxwellian limits. Solutions of the wave- 
dispersion relation show that these modes may be unstable for some plasma parameters. 
Both convective and absolute instabilities were found near the cyclotron frequency whereas 
only convective instabilities were found at the cyclotron harmonics. In addition, it was 
found that there could be present an absolute instability near the hybrid frequency and 
convective instabilities at frequencies between the cyclotron harmonics. It is well known 
that cyclotron harmonic waves may be unstable in the collisionless limit for some of the 
electron distribution functions considered in this section. However, all cases investigated 
here are stable or evanescent in the collisionless limit; therefore, the electron-neutral 
collision process is indeed responsible for the predicted instabilities. 

THEORY OF WAVES PROPAGATING PARALLEL 
TO THE MAGNETIC FIELD 

In this section waves propagating along the direction of the applied magnetic field 
are considered for both isotropic and anisotropic collision models. The general case of 
anisotropic collisions is considered initially and the first-order distribution function is 
written in terms of a set of inhomogeneous partial difference equations. After special- 
izing to waves propagating parallel to the magnetic field, these equations are solved and 
the conductivity tensor and dispersion relations are obtained. The isotropic collision 
model is then considered and the conductivity tensor and dispersion relations are again 
derived. It is shown that in the limit of isotropic collisions, the dispersion relations 
found from the anisotropic collision frequency model are equivalent to those found when 
collisions are assumed to be isotropic in the beginning. Some of the analytic properties 
of the dispersion relations are discussed and comparison is made with other collision 
models. 


First-Order Distribution Function With Anisotropic Collisions 

In the Lorentz gas model it is assumed that the ions and neutrals are infinitely 
heavy and that the plasma is weakly ionized so that electron-neutral collisions are much 
more frequent than either electron-electron or electron-ion collisions. The Boltzmann 
equation for electrons can be written as 

- + V ~ - Me + v X §) • ^ = N|v I C ct(v,x) [f(v ') - f(v to (79) 

at a r m av J u - 1 
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where x is the angle between v and v\ The number density of the neutrals is N, 
and a(v,x) is the differential scattering cross section describing the electron-neutral 
collision process. This equation is linearized by making the ansatz 


f(v) = f 0 (v) + f x (v) 

E = 0 + Ej \ 

B = B 0 + Bf v 


(80) 


The zero-order equation 

9f 

(v x W c ) • = N|v | y ct(v,x) [f 0 (v«> - f 0 (v)jdf2’ 

and the first-order equation 


9f t af af, Q x af n 

— + v • — - (v xS.) • -J- = — (e, + v x b. ) • -J. 

at f)r ' c ' av ni\ i -*■' av 


n|v| j* a(v,x) [f^v’) - f 1 (v)|dJ2' 


(si; 


(82 


where a> c = eB^^/m are then obtained. The solution of the zero-order equation 
(appendix A) is any isotropic distribution function fg(v) and hence the term 
(v X B^ • df^/dv in the first-order equation is identically zero. As done previously, 
the Fourier transform in space and the Laplace transform, in time, corresponding to the 
wave representation exp(iu>t - ik • ?) are taken to obtain 


*af - _ 9f n 

i(u> - v" • k)f. - (v x o> r ) • = — E n • — - 

IV cj m 1 


+ N|v | j 1 ct(v,x) I?! (?’) - fj^dn’ 


Expanding f ^(w,k ,v) in spherical harmonics in velocity space yields 

°° l 


1=0 m=-f 


(83) 


(84) 
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, 1/2 


where Y, (0,<p) = ~ ”^’1 P 7 m (cos 0)e im<5 °. Substituting this relation into 

L,m ^ 47j-(£ + m). J 4 

equation (83) yields 


I i (“ - 7 ' * + * §■ ®i • - N|v| J ^ m Yj <% 

Z,m dv l,m 


m 


(85) 


where 


Q i,m ■ y ' 7,v '*' 




df!' 


In appendix J it is shown that for all m, 

Q^ m =( ^i = 2n £ CT ( V >X) [* " P^cos x)] sin x dx 

where P^(cos x) is the Legendre function of degree l. Then by following Allis 
(ref. 42) and Derfler (ref. 14), the transfer collision frequencies are defined as 

^(v) = N|v |Qj 

and equation (85) then yields 

V — — 9f n 

) i(u> - iv 7 - v • k + mwjf, Y, = — E., • 

Lu \ t wZ,mZ,m ml 9v 

Z,m 


(86) 


(87) 


( 88 ) 


(89) 


Multiplying this equation by Y* „(6,<P) and performing an integration over the solid 

r,s 


angle dO = sin 9 d 0 d cp, an inhomogeneous difference equation is obtained 

i(o> - iv r + sw c )f r>s - i £ i l m y v . k\ m (0,<p) y; > 8 (0,^) df2 = T r p (90) 

Z,m 


where 


l/2 ”| 

T r,s = f f i(x) 6 r,lj^s,l( -E x + iE y) + 6 s,-l( E x + iE y) + 6 s,0'^ YE z] <«> 

The integral can be easily evaluated for general k to give a partial difference equation 

for f_ but instead it will be specialized to the case of k parallel to the magnetic 
r,s 
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field. For this case, v • k = vk cos 9 and the integral becomes 

V, y cos . y* «. -^f | rts ^lV S + 1) ] 

1/2 


1/2 


5. 


Z,r+1 


~(r + s)(r - s )~| 1/2 6 

2r - 1 J V-l 


From equation (90) the difference equation 


i/o) - H'r + a« c )f r s Rr + s 4 JHr - s + 1)' 

t r c) r ’ s ^ r+1 ’ s L 2r + 3 J 

. f Rr + s)(r - s)f 
r -!, s L 2r - 1 


1/2 


( 92 ) 


(93) 


is obtained. Note that this equation is an ordinary difference equation in the index r at 
constant s and is readily solvable by use of standard techniques. This equation is 
solved later in this section to obtain the conductivity tensor and dispersion relation. 

For the case k 1 B, equation (90) would lead to a partial difference equation, that 
is, one in which both the indices change. This type of equation is not readily solvable by 
using standard techniques and is not considered in this report. Note, however, that the 
isotropic collision model for k 1 B was treated earlier. 


Conductivity Tensor and Dispersion Relations for 
k || B With Anisotropic Collisions 

The conductivity tensor fof the case of propagation parallel to the magnetic field 
will now be sought. Equation (84) is substituted into the current equation 


j = a • E = -en 


y vf-Rv^d*" 


(94) 


to obtain 


v3 R( f i,i- f i,-i) + V( f i,i +f i,-i) - 5 c' / 2 f i,o] tiv 


(95) 


Note that only f ^ _j, f ^ q, and f ^ ^ are needed to find the conductivity tensor d, and 

with this in mind the solution of equation (93) is expressed in terms of a continued fraction 
as 
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(96) 


f l,s g s T l,s 


where 


g S = 


3/ikv 4 - s 2 9 - s 2 16 - s 2 
+ 3a T,s “ 5X 2,s " 7X 3,s " 9A -4,s 


n 2 - s 2 


(2n + 1)A. 


n,s 


(97) 


and 



(s = ±1) 
(s = 0) 


n,s 


co - IV + SO) 

n c 

kv 


(98) 


(99) 


s = 0, +1, -1 

The continued fraction notation used here is defined in references 43 and 44. By combining 
equations (95) and (96), the conductivity tensor may be written as 


'll 


a ~ 


-a. 


12 


12 


’ll 


0 

0 

T 33 y 


( 100 ) 


where 


CT 11 “ 


°12 = 


°33 “ 


• c o"pfI 0 v 3 f o( g i + g -i) dv 
• e o"pfl 0 vS, o i (*-i - g l) dV 
■**•$7 S'**** 


( 101 ) 

( 102 ) 

(103) 
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From appendix C, with k = keg in equation (C17), the dispersion function is 


D(w,k) = 


i + ( a n / iCt , e o ) ~ (<*/“>)' 

12/ iwe 0 

0 


-a. 


a iih e o 

1 * ( II ll / i " < o ) ■ (ck ' /, -'-' )2 
0 


0 

0 


1 + 


( < T 33 / 1 " e o ) 


= 0 


(104) 


Expansion of the determinant yields two uncoupled dispersion relations 


Uoo 

Dp(w,k) H 1 +T-^-= 0 
E' 1CUC q 


(105) 


**<".*> ■ i - (f ) 


2 ( J 11 * ig 12 _ „ 
iu)e n 


(106) 


By using equation (103) for o^g, equation (105) becomes 


D F (w,k) = 1 + -£ C v 3 f’g n dv = 0 


3 u> Jr 


0 


f oSo 


(107) 


This dispersion relation represents a purely longitudinal electrostatic wave that is not 
affected by the magnetic field. This same dispersion relation was found by Derfler 
(ref. 14) for an unmagnetized plasma and is not considered further here. Substituting 
for Oji and * n equation (106) yields 


/ „ b .v2 2iru* r»°° 

D ± (u>,k) = 1 - (-1) + 3 o g ± (v) v2f 0 (v)dv = 0 


(108) 


where g is the continued fraction 


g * W = 3\ 


n 2 - 1 


1,±1 “ 5 X 2,±1 


- (2n + l)X n>±1 - 


(109) 


This dispersion relation describes transverse electromagnetic waves propagating parallel 
to the applied magnetic field, that is, left- and right-hand circularly polarized waves. 
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The dispersion relations derived in this section (eqs. (107) and (108)) exactly 
describe the dispersion of waves propagating along the magnetic field in a Lorentz gas 
plasma. If an anisotropic electron-neutral collision cross section cr(v,x) is given, one 
can obtain the transfer collision frequencies from equations (87) and (88) and then, in 
principle, solve the dispersion relations. Often, in practice, only the velocity dependence 
of the collision cross section is available and in this case the isotropic collision model is 
of more interest. This case is considered in the next section. 


Conductivity Tensor and Dispersion Relations for 
k | B With Isotropic Collisions 

In this section it is assumed from the beginning that the electron-neutral collision 
process is isotropic. This procedure was followed earlier for general k where the 
calculation started with the Boltzmann equation given by equation (14) and a conductivity 
tensor a given by equation (34) was found. In appendix K this a is specialized to 
propagation parallel to the magnetic field to obtain 



ct 12 

a ll 

0 



where 


a n = - 1C V 


r f o v i 

0 J ~v~~ 


U) 


kVr 


H 3 
a v 


c o r 


■V “ ' 3 ) 


fw - kv. 


'12 


p 0 J V 


w. 


Cl) 




d 3 v 


CT 33 = ^ 2 — 





kv 


(HO) 


(111) 


( 112 ) 


(113) 


= co - it-(v) 
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After writing d^v = sin 9 d0 d (p dv, the angular integrals in equations (111) and (112) 
can be evaluated exactly but it is easier here to leave them in their present form for now. 
Using equation (105) and Ogg from equation (113) yields 


D E (w,k) = 1-47 t 



(114) 


where 


X 


LO v 

kv 


This equation describes a purely longitudinal electrostatic wave as does equation (107) 
and Derfler (ref. 14) has shown that in the limit of isotropic collisions, equation (107) 
reduces to equation (114). 

From equations (106), (111), and (112), the following dispersion relation for trans- 
verse electromagnetic waves is obtained: 




f 0 v l/ v 


Up - kv 3 ± w c 


d v = 0 


(115) 


9 9 

Then since d v = v^ sin 9 d0 d cp dv, the integrals over the angles can be evaluated to 
obtain 


D ± ( W ,k) = 1 - (f ) 2 + v2f' \x ± + (l - X^coth- 1 


dv = 0 


(116) 


where 


X ± = 


CO ± (jO 

w V c 
kv 


(117) 


In appendix L it is shown that equation (108) is identical to equation (116) if isotropic col- 
lisions are assumed in the former. In appendix M equation (116) is derived in a much 
simpler but less rigorous manner. The dispersion relation expressed by equation (116) 
is discussed further in the next section and is solved subsequently. 
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Discussion of the Transverse Wave Dispersion 
Relation for k II B 

If one goes through the derivation of the dispersion relation (eq. (116)) carefully, it 
is apparent that the integral term in equation (14) has no contribution as was shown in 
reference 45. This condition is true only for an isotropic collision process and for 
transverse waves propagating parallel to the magnetic field. Hence, the collision model 


(SLr - ,(5r) ] 


(118) 


will give the identical dispersion relation. If there is no velocity dependence of the col- 
lision frequency in equation (118), then the constant collision frequency model of Allis 
is obtained. Hence, the dispersion relation for constant collision frequency is obtained 
merely by setting v equal to a constant in equation (116). This model is often used 
because of its relative simplicity both in derivation and solution. In the case of a 
Maxwellian electron distribution function, that is, equation (75), it is found in appendix N 
that this model yields the dispersion relation 


D ± (o>,k) = 1 




wkvQ 


Z 



, ±U) c\ 

kv 0 / 


= 0 


(119) 


where Z^Aq^ can be written in terms of the complex conjugate of the Fried function 
(ref. 35) (plasma dispersion function) by 


Z 




( 120 ) 


This dispersion relation has been previously obtained by Scarf (ref. 46), Lee (ref. 47), 
and others. 

In order to show how the dispersion relation given by equation (108) or (116) deter- 
mines the propagation characteristics of waves, it may be assumed that there is present 
a localized driving force given by 



( 121 ) 


Substituting this equation and equation (110) into Maxwell's equations and performing 
appropriate manipulations yields 


E ± (o>,k) 


ij ± Kk) 

Ct)e 0 D T (o),k) 


(122) 
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where 


E ± = E x ± iE y 


and 


j ± i x ± ijy 

Hence, it is seen that D (w,k) and D + (a>,k) represent right- and left-hand polarized 
waves, respectively. 

When solving the dispersion relation, it is important to keep in mind that equa- 
tion (116) is defined for k real and o> on a Laplace integral path. For this case the 
real v-axis is the proper contour as indicated. If solutions for k(u> real) are desired, 
it is necessary to use the appropriate analytic continuation of D(o>,k) which is found by 
deforming the contour in complex velocity space so that no singularities cross the con- 
tour while u> is moved from the Laplace integral path (LIP) to real frequencies. 

Note that the symmetry relation 

D + (w,k) = D*(-u>*,-k*) (123) 

is satisfied, and hence it is sufficient to consider only one of the dispersion relations: 

D (w,k) =0 or D + (oj,k) = 0. Also note that 

D(u>,k) = D(w,-k) (124) 

thus, k = 0 is a saddlepoint in the k plane, and hence the roots of D(w,0) = 0 will be 
branch points in the u> plane. The location of these branch points in the frequency 
plane is of prime importance in determining the stability character of the system 
(refs. 23 to 26). Details are given in the next section where solutions of the dispersion 
relation (eq. (116)) are found. 

SOLUTIONS OF THE DISPERSION RELATION FOR TRANSVERSE 
WAVES PROPAGATING PARALLEL TO THE MAGNETIC FIELD 

In this section the dispersion relation for propagation parallel to the magnetic field 
is solved for several different electron distribution functions. A new collisional mode is 
found to be unstable for some ranges of the plasma parameters. The wave considered is 
the transverse -electric, circularly polarized wave known in ionospheric work as the 
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"whistler" wave. In solid-state plasmas it is sometimes called a "helicon" wave. Both 
convective and absolute instabilities are found when the electron distribution function is 
monoenergetic. A Lorentz gas with nitrogen neutrals is then considered and instabilities 
are again found when a peaked electron distribution function is considered. A Maxwellian 
electron distribution function is then assumed, and it is found that the wave dispersion 
predicted by the isotropic collision model (ref. 48) is very close to that predicted by the 
constant collision frequency model of Allis. An electron distribution function consisting 
of both peaked and Maxwellian components is briefly considered, and the stabilizing 
influence of the Maxwellian population is demonstrated. 


Monoenergetic Electron Distribution 

For the initial analysis of the dispersion relation (eq. (116)) it is assumed that the 
zero-order electron distribution is monoenergetic and solutions of the dispersion rela- 
tion are obtained. Substituting equation (66) into equation (116) where the negative sign is 
chosen so that the dispersion relation represents right-hand polarized waves yields 


D(w,k) = 1 


m 


w p 

iliy 

wkvQ 

kv 0 


(l - X coth - lx) - 


coth *X 


= 0 


where 


(125) 


h = 



v=v 


0 


CO, 


CO 


X = 


kv r 



This dispersion relation was also derived and solved by Suzuki (ref. 15) but here a some- 
what more detailed analysis of the stability of the waves is given. 

Note that coth~l\ is a multivalued function of X; therefore, equation (125) has an 
infinite set of roots. However, later it will be shown that only those roots on the principal 
sheet are of physical interest. As Vq goes to zero, equation (125) approaches the 
familiar cold plasma dispersion relation 


D(co,k) 


1 - 



CO 


2 

P 


CO f CO 



= 0 


(126) 


consisting of the upper and lower electromagnetic fast-wave branches and the slow-wave 
whistler mode. 
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The dispersion relation given by equation (125) has been solved under conditions 
approximating those in the ionosphere at an altitude of 80 km, that is, u> p = 2u> c and 
v = 0.06u> c . Nitrogen exhibits a resonance scattering (ref. 39) such that h ~ 5 at 

mvQ^2e =1.5 eV. By using these values, the results of figure 14 that show the real part 
of the complex frequency as a function of k real are obtained. The upper, lower, and 
whistler branches (U, L, and Whl, respectively) deviate only slightly from their cold 
plasma values. The Wh2 mode is a new collisional mode not predicted by cold plasma 
theory and first recognized by Suzuki (ref. 15). It has a very small group velocity and 
hence represents a plasma resonance. When either h or v p approaches zero, this 
mode has the solution w = oo c .+ iv for all k. The major difference between the solu- 
tion in figure 14 and the cold plasma solution is the mapping of the roots, Whl and Wh2, 



kc/o; c 

Figure 14.- Brillouin diagram, Re(cu) plotted against k real. 
Dotted line represents higher order sheet. 
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onto higher order Riemann sheets as k becomes sufficiently large. This mapping is 
shown for the mode Whl in figure 15 where the real k axis is mapped into the com- 
plex (jo plane by the dispersion relation with h = 0. In order to understand the 
behavior of this root, note that coth“^X has a branch cut connecting the points at 
X = ±1. If as k varies, the co root moves in such a manner that X crosses the real 
axis between +1 and -1, the appropriate analytic continuation of coth*“^A. must be used. 
By writing co in terms of X, 

co = co c + iv + Xkvg (127) 

and taking the branch cut to be on the real axis in the X plane, it is seen that in the 
co plane the cut for k real goes from Re(cu) = a> c - kvQ to Re(o>) = o> c + kv Q along 
the line Im(w) = v. Hence, whenever the c o root crosses this line, X must go onto 
another Riemann sheet. For the case shown in figure 15, it was found that the root 
oscillates with extremely small excursions about the branch line Im(w) = v and thereby 
goes onto successively higher order Riemann sheets. However, once the root leaves the 
principal sheet, it no longer needs to be considered explicitly, since in constructing a 
Green’s function, this root would be taken into account by an integration around a branch 
cut. This phenomenon was explained in more detail by Derfler (ref. 14) for the case of 
electrostatic waves in an isotropic Lorentz gas. 



Figure 15 Locus of real wavenumber k in the frequency plane for 
mode Whl. Dotted line represents higher order sheets. 
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In figure 16, results are presented for the case h = ±3 and h = ±5 where the 
real k axis is mapped into the complex frequency plane for both the Whl and Wh2 modes. 



Figure 16 .- Locus of real wavenumber k in the complex frequency plane 
for modes Whl and Wh2. Dotted line represents higher order sheets. 

In figure 16(a) the roots never dip into the lower half plane; thus, it is concluded that no 
instabilities exist for that case. In figure 16(b) it is seen that for h = 5, the mode Wh2 
has a branch point with Re(o>) ~-o> c in the lower half plane and indicates an absolute 
instability, whereas for h = -5, the mode Whl dips into the negative half plane and indi- 
cates a convective instability (refs. 23 to 26). A similar mapping for the fast electro- 
magnetic modes in figure 17 shows that both branches become absolutely unstable for 
h < -3. 

Negative values of h less than -3 are seldom realized in nature, and therefore, 
the new collisional mode Wh2 is considered more important for this analysis as it goes 
unstable for values of h 5 3. Since a value of h = 5 is characteristic of the resonance 
scattering peak near = 1.5 eV in nitrogen, h = 5 was chosen to investigate 

the growth rate Im(cu) of the Wh2 mode at k = 0. The growth rate is given by the 
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Figure 17*- Locus of real wavenumber k in the complex frequency 
plane for modes U and L. 

solution of equation (125) with k = 0, that is, by 

w - ^(W3) (co p - u> c y 2 - - u^)- 1 = 0 (128) 

In figure 18, the growth rate Im[u>(k = 0)[] is shown as a function of v with u>p as a 
parameter. Instabilities with an e-folding time of 1 cyclotron period are easily obtained 
with the parameters Wp > u> c and v > 0.2co c . Such values can be obtained in laboratory 
plasmas in the range of 0.1 to 10 mm Hg (13.33 to 1333 N/m^). The lower left-hand 
corner of the diagram is typical of the conditions in the lower ionosphere in which devia- 
tions from a Maxwellian equilibrium are known to occur (refs. 49 and 50). It is therefore 
possible that collisional instabilities of this type arise in nature. 

Equation (128) can be rearranged to yield a third-order polynomial in oo. The 
exact locations of the roots will depend on the plasma parameters chosen but for a large 
variety of cases it is found that the real parts on the roots were near u> c and 

w c/ 2 ± H/ 4 + w 2 ) ^ . If either v or h goes to zero, the order of equation (128) 
is reduced by 1. The root with real part near u> c vanishes, and hence it is concluded 
that this mode was introduced by collisions. 
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Figure 18.- Growth rate Im(<x>) of Wh2 mode at infinite wavelength as a 
function of collision frequency v. Re(co) =* 0 ^. 


Peaked Electron Distribution 

Equation (116) was solved for a monoenergetic electron distribution function and 
instabilities were found for certain ranges of the slope parameter h. In particular, it 
was found that for h > 3, there was an unstable mode (Wh2) with oj ~ u> c . For h < -3 
three unstable modes, a whistler mode (Whl) and the two fast modes (U and L), were 
found. Hence, one concludes that collisional instabilities of this type will be found in a 
plasma if there are present (1) a background of neutrals that has |h| > 3 over a region 
or regions in velocity space and (2) a sharply peaked electron distribution function cen- 
tered in one of these regions. As was done previously, the experimental electron-neutral 
collision cross section for nitrogen (eq. (72)) is used. To meet the second criterion, 
again the peaked distribution function described earlier and given by equation (73) was 
used. 

A computer program was written to solve the dispersion relation (eq. (116)). The 
velocity integral was performed numerically by using a Gauss quadrature rule. (See 
appendix G). The integrand has a branch cut between the points \ = 1 and X = -1, and 
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upon mapping this line into velocity space, six distinct branch cuts are found. The 
proper contour is found by requiring that the real v-axis from zero to infinity be a valid 
contour when k is real and when w is on a Laplace integral path (LIP). Generally, 
the real v-axis from zero to infinity was found to be a valid contour for other u> and k, 
but occasionally it was necessary to deform the contour as described in appendix I. 

The real k-axis was mapped into the complex frequency plane (fig. 19) to determine 
the locations and types of instabilities. As expected, the new collisional mode was 
unstable when s was sufficiently large, and when Vq was located on the positive slope 
of the collision frequency curve. No instabilities were found when Vq was located on 
the negative slope of the collision frequency curve, as numerical difficulties were 
encountered before the half -width of the distribution function could be made small enough. 
The instability of the mode with Re(cu) ~ co c was investigated more fully, and the curves 
in figure 20 giving the boundaries in the s,Vq plane between the stable, convectively 
unstable, and absolutely unstable cases were obtained. For the points A and B in fig- 
ure 20, the real k axis maps into the complex frequency plane as shown in figure 21. 

In figure 22 the electron distribution function is shown plotted for the points C and D in 
figure 20. 

The value of the collision frequency chosen for these figures, that is, = 5.26o) c , 
represents conditions in the lower ionosphere. Note that is the collision frequency 



Figure 19.- Real k axis mapped into the complex frequency plane 
for peaked electron distribution. Plot shows absolutely 
unstable waves. 
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Figure 22.- Normalized peaked distribution function for points on 
boundaries of regions in figure 20. Points C and D. 

at infinite electron velocity. If the collision frequency is averaged over a realistic elec- 
tron distribution function, the average collision frequency will be on the order of 2u> c . 
This value is perhaps on the high side of the range of values found in the ionosphere but 
lower values of will only change the strength of the instability and will not, for 
example, cause the case presented in figure 19 to become stable. It should be noted, 
however, that in the collisionless limit the collisional mode at w ~ u> c in figure 19 
vanishes, and the plasma is stable. Since the plasma is also stable in the collisional 
limit, there is some intermediate v x (and hence a neutral particle density) at which the 
collisional mode is most unstable. 

For either stable or convectively unstable waves, the mapping of real o> into the 
wavenumber plane is of interest as this case is the steady- state situation often investi- 
gated in the laboratory. In figures 23 and 24, Brillouin diagrams that result from such 
mappings are shown. For comparison, the prediction of the constant collision frequency 
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Figure 23.- Brillouin diagram giving Re(k) ard Im(k) as a function 
of cd real for a convecti-rely unstable case. 

model of Allis obtained by using equation (10) for the collision term is also shown. The 
value of the collision frequency used in this model was the average over the distribution 
function, that is, equation (77). The exact theory predicts a convective instability for the 
case presented in figure 23 since Im(k) goes positive near u> = u> c . 

Figure 24(a) is a similar case that is stable because of the slightly higher Vq. 
Figures 24(b) and 24(c) are similar to figure 24(a) except that the collision frequencies 
are, respectively, one-fifth and one-tenth of those in figure 24(a). In figure 24(c) the 
collisionless limit is also shown and from the sequence figure 24(a) to figure 24(c) it is 
possible to see how the curves are approaching the proper collisionless limit. 

Comparison of the Allis model and the exact theory shows that the discrepancies 
are rather large. Particularly near u> = o> c , the exact theory predicts decreased 
damping and, in some cases, instabilities. Hence, it is seen that for non -Maxwellian 
electron distribution functions, it is necessary to treat collisions rigorously to obtain the 
correct wave dispersion. 
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(a) v oo = 5 . 26 a> c ; <v> = 2.0^,. 

Figure 24.- Brillouin diagram giving Re(k) and Im(k) as a function 
of CD real for a stable case . 








Maxwellian Electron Distribution 


A Maxwellian electron distribution (eq. (75)) was then assumed and equation (116) 
was solved by use of the same collision frequency function. Since no instabilities were 
expected, the real frequency axis was immediately mapped into the wavenumber plane to 
obtain the Brillouin diagrams shown in figures 25(a) and 25(b). Again, the predictions of 
the constant-collision-frequency Allis model are shown for comparison purposes. The 
average of the velocity-dependent collision frequency as defined by equation (77) was 
used in this model. 

The agreement between these two models is clearly good. Several other cases 
were calculated and the agreement between the two models was similar to the cases pre- 
sented here. Hence, it would seem that if a Maxwellian electron plasma is being con- 
sidered, the dispersion relation given by equation (119) should be satisfactory for most 
purposes. One could then use the very fast existing computer program for the plasma 
dispersion function developed by Derfler and Simonen (ref. 51). 



Im ( kc/cu c ) Re ( kc /cu c ) 

(a) v Q = l.iH\/eV; v M = 5 . 260 ^.; <v> = 1.42o^. 

Figure 25.- Brillouin diagram giving Re(k) and Im(k) as a function 
of cd real for a Maxwellian electron distribution . 
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Sum of Maxwellian and Peaked Electron Distributions 

Since it was found that a peaked electron distribution may be unstable, one naturally 
wonders how great would be the stabilizing influence of a Maxwellian electron population 
when added to a peaked distribution. To answer this question, equation (78) was again 
used as the electron distribution function. For the peaked component, parameters were 
used that were known to lead to an absolute instability, that is, those used in figure 19. 

The stability of the waves at various values of a and v^ was then investigated. The 
results are shown plotted in figure 26. From this figure it can be seen that the stabilizing 
influence of the Maxwellian population is critically dependent on the temperature of the 
Maxwellian. Since absolute instabilities occur at zero group velocity, that is, 9u>/ 9k = 0, 
it is not surprising to find that low-temperature Maxwellians have more of a stabilizing 
influence because of the larger number of low-velocity electrons. Comparison of fig- 
ures 26 and 12 shows that the Maxwellian population required for stabilization is much 
less for transverse waves propagating parallel to the magnetic field than for longitudinal 
waves traveling perpendicular to the magnetic field. In figure 27 the electron distribution 
function for a case on each boundary in figure 26 (points A and B) has been plotted. Note 
that both cases are very much dominated by the peaked population of electrons. 
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MAXWELLIAN POPULATION IN PERCENT 



THERMAL SPEED v, OF MAXWELLIAN IN ,/eV 

Figure 26.- Stability plot for peaked plus Maxwellian distribution 
function. Mode with Re(a>) ~ cd c . 
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Figure 27.- Normalized distribution function consisting of peaked and 
Maxwellian components (eq. ( 78 )). Points A and B in figure 26. 

Synopsis of Numerical Results 

Briefly, it has been found that a rigorous treatment of electron-neutral collisions in 
a non-Maxwellian Lorentz magnetoplasma predicts a new mode as well as the usual modi- 
fied collisionless modes. This new mode vanishes both in the Maxwellian and collision- 
less limits and may cause the plasma to be unstable under some conditions. This insta- 
bility may be either convective or absolute and occurs at frequencies near the cyclotron 
frequency. The modified collisionless modes were also found to be unstable for some 
plasma parameters. It is known (ref. 47) that in the collisionless limit the transverse 
electromagnetic waves considered in this section are stable; therefore, the instabilities 
investigated here are a result of the electron-neutral collision process. 
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DISCUSSION 


The purpose of this report was to determine rigorously the effects of collisions on 
electron waves in a weakly ionized magnetoplasma. The waves were assumed to propa- 
gate in an electron plasma with infinitely heavy neutrals and ions. By using this model 
(usually called the Lorentz gas model), the linearized Boltzmann equation was solved 
with the collision integral and the conductivity tensor was derived. The dispersion rela- 
tions for waves propagating both perpendicular and parallel to the magnetic field were 
then derived from Maxwell's equations and solved for a number of cases of general inter- 
est. It was found that the plasma supported a new class of collisional modes if (1) the 
collision frequency was velocity dependent, and if (2) the electron distribution function 
had a region in velocity space where 9f q/ 9v > 0. These new modes could be unstable if 
a large proportion (> 70 percent) of the electrons were located in a region in velocity 
space where the collision frequency increased (or decreased) rapidly. 

For waves propagating perpendicular to the magnetic field, that is, k 1 B, the col- 
lision model assumed allows for an isotropic velocity -dependent electron-neutral colli- 
sion cross section, and two uncoupled dispersion relations were derived. One described 
a modification of waves known in cold plasma theory as ordinary waves since they are 
unaffected by the magnetic field in that limit. The other described extraordinary waves 
and has slow wave solutions known variously as cyclotron, Bernstein, electrostatic, or 
longitudinal waves. In the collisionless limit, this dispersion relation yields an equation 
equivalent to that of Bernstein (ref. 29), whereas in the limit of zero magnetic field 
Drummond's results (ref. 13) are recovered. When both the collisions and the magnetic 
field vanish, Landau’s dispersion relation is obtained (refs. 32 and 33). 

The collisional dispersion relation for k 1 S was solved for several electron dis- 
tribution functions. One series of solutions were damped modes that were modifications 
of the undamped collisionless electrostatic modes. Under some conditions they could be 
driven unstable at frequencies near zero and the hybrid frequency. In addition, a series 
of new collisional waves were found with real frequencies near zero and near each cyclo- 
tron harmonic. Under the proper conditions, either absolute or convective instabilities 
were possible at the cyclotron frequency and convective instabilities were possible at the 
harmonics. 

Waves propagating parallel to the magnetic field were also considered, and a colli- 
sion model allowing a velocity and an angle-dependent electron-neutral collision cross 
section led to two uncoupled dispersion relations. One described an electrostatic or 
longitudinal wave that was not affected by the magnetic field. This same equation was 
derived by Derfler (ref. 14) for an unmagnetized plasma. The second dispersion relation 
described a transverse wave, one mode of which is commonly called a whistler wave in 
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ionospheric work. This dispersion relation was solved and a new collisional mode was 
found near the electron cyclotron frequency. It was found that this mode could be either 
absolutely or convectively unstable under some conditions. 

The isotropic electron distribution functions used when solving the dispersion rela- 
tions were of varying physical interest. The monoenergetic distribution function was of 
interest primarily because it gave relatively simple equations to solve. Furthermore, 
the solutions were very helpful in showing the type of results to expect when more real- 
istic distributions were used. The peaked, Maxwellian, and Maxwellian with a bump on 
the tail were all of much more physcial interest since these distributions are quite pos- 
sible in various types of plasmas. 

To actually observe and measure the wave instabilities that this theory predicts, 
some means has to be devised to set up a known electron distribution function with a peak 
of the proper width and location. The plasma could then be excited at some frequency 
and the complex wavenumber measured. It would not be necessary, however, to actually 
observe an instability to show the validity of this theory, since it has been shown that for 
non-Maxwellian plasmas, the predicted wave dispersal may differ greatly from that pre- 
dicted by simpler collision models. For example, experiments that have been previously 
done in the Vlasov limit might be redone in a regime where collisional effects are impor- 
tant. For the electrostatic waves with k JL B, one should look for decreased damping at 
frequencies just below the cyclotron harmonics. For the transverse waves with k || B, 
the decreased damping should occur precisely at the cyclotron frequency. If these effects 
are observed, it should be possible to assume some model for the electron distribution 
function, for example, equation (78), and vary the parameters to obtain the best fit to the 
experimental wave dispersion data. This procedure was followed by Fields et al. (ref. 3) 
using an equivalent temperature approach. 

Most experimenters to date have merely observed the radiation from the plasma in 
an attempt to observe the collisional instabilities. Tanaka and Takayama (ref. 5) and 
Oddou (ref. 4) and others have observed anomalous emission at the cyclotron frequency 
in plasmas containing the three noble gases: argon, xenon, and krypton. These gases 

are generally called the Ramsauer gases, and they have a region where the electron col- 
lision frequency curve is very steep. No anomalous emission was observed from plasmas 
containing the other two noble gases: helium and neon. These gases have relatively flat 
collision frequency curves, and the theory given in this report would not predict any col- 
lisional instabilities. The emission from the Ramsauer gases was called anomalous 
since it had an intensity many orders of magnitude larger than could otherwise be 
explained. These experimenters observed also that the emission occurred in bursts that 
had a direct relationship to an oscillation of much lower frequency. It seems likely that 
the low-frequency oscillation was setting up the appropriate electron distribution function 



necessary for the plasma to be unstable at the cyclotron frequency. It should be noted 
that the radiation was observed primarily at the electron cyclotron frequency. The 
theory presented here predicts an absolute instability at the cyclotron frequency but not 
at the higher harmonics. Hence, much stronger instabilities (and hence more radiation) 
would be expected at the cyclotron frequency than at the harmonics. 

Wachtel and Hirshfield (ref. 8) have also performed an experiment demonstrating 
collisional instabilities. They introduced monoenergetic electrons into a microwave 
cavity containing xenon. They observed negative absorption at the cyclotron frequency 
when the speed of the electrons was in a region in velocity space where the collision 
frequency had a large positive slope. By negative absorption they mean that the electrons 
were giving up energy to the waves in the cavity, and hence the waves were unstable. 

The collision-induced instabilities investigated quantitatively in this report can be 
understood qualitatively by considering a very simple model. Take a wave of frequency 
w ~ propagating in a magnetoplasma consisting of monoenergetic electrons. For the 
case of the transverse electromagnetic waves, right circularly polarized waves with 
k || B are assumed, and for the electrostatic waves one has iT 1 B with E II k. Assume 
that there are N electrons in the plasma. A little thought will show that in the begin- 
ning for either direction of propagation, N/2 electrons will be gaining energy from the 
wave and N/2 electrons will deliver energy to the wave. After a given length of time, 
At, let Cq be the number of electrons that had a collision with a neutral while gaining 
energy. Of those electrons, Cq/ 2 will be gaining (or losing) energy after the collision 
since the collision will tend to randomize momenta. Likewise, let Cq be the number 
of electrons that have a collision while losing energy. After time At, the number of 
electrons that are gaining energy is 






and the number of electrons that are losing energy is 





(129) 


(130) 


If more electrons are losing energy, then the wave must be gaining energy and hence is 
unstable. Hence, assume that Nq > Nq; equations (129) and (130) then yield Cq > Cq. 
But the number of electrons that collide in a given time is proportional to the collision 
frequency, thus, viy + Av) > v(y - Av). Hence, the collision frequency must have a 
positive slope for the plasma to support an unstable wave at the cyclotron frequency. 
Earlier, it was found that these wave instabilities were indeed present only if the slope 
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parameter was positive and large enough. Because of the extreme simplicity of the 

model considered here, it is not surprising that the agreement is qualitative in nature 

■6 

and not quantitative. It is also interesting to note that this simple argument predicts 
instabilities for cv = o? c but not for w = nu> c (n = 2,3,4,. . .). It may be recalled that 
solutions of the dispersion found earlier sometimes predicted absolute instabilities at the 
cyclotron frequency but never at any of the harmonics. 

The theory given in this report appears to be the most accurate and extensive to 
date concerning the effects of electron-neutral collisions on high frequency waves in a 
magnetoplasma. The model assumed is that of a plasma of electrons colliding only with 
infinitely heavy neutrals, and since all ion effects are ignored, this theory is valid for 
frequencies somewhat above the ion plasma and ion cyclotron frequencies. By using the 
techniques described in this report, it is possible to solve the dispersion relations for 
essentially any electron distribution function, and by working backward one might be able 
to calculate the electron distribution function or the collision frequency function when 
given the experimentally measured wave dispersion. The equations derived and the 
techniques described in this report should be particularly applicable to wave propagation 
through and radiation from plasmas containing especially the Ramsauer gases (Ar, Xe, 
and Kr) since the electron-neutral collision cross sections in these gases are such that 
large collisional effects on waves are both expected and observed. 

CONCLUDING REMARKS 

The Boltzmann equation with collision integral was solved and the wave dispersion 
relations for electron waves propagating both parallel and perpendicular to the magnetic 
field were obtained. A Lorentz gas magnetoplasma was the model assumed. This model 
requires (1) that the ions are infinitely heavy and (2) that the electrons are colliding only 
with infinitely heavy neutrals. Numerical solutions to the wave dispersion relations were 
obtained, and for both directions propagation the electron waves could be unstable if a 
large fraction of the electrons had a velocity in a region where the electron-neutral 
collision frequency increases (or decreases) rapidly enough. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., May 4, 1971. 
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APPENDIX A 


SOLUTION OF ZERO-ORDER BOLTZMANN EQUATION 

In this appendix it is shown that any isotropic distribution f g ( | v | ) is a solution of 
the zero-order Boltzmann equation (eq. (16)) which in spherical coordinates becomes 


w. 


3f 0 (v) 

d<p 


= + ^ y 


(Al) 


where Bg is taken equal to BgCig. Let us write fg(v) as a spherical harmonic 
expansion in velocity space 


00 ^ 

l - 0 m=-Z 


L 

f 0 W = I 2 


(A2) 


Substituting equation (A2) into equation (Al) yields 




7,m v 


(A3) 


For a nonzero v, this relation is satisfied only by f^ m = 0 for all l and m except 
l = m = 0. Hence, from equation (A2), 

f 0 (r) =f o,o (v) ( A4 ) 

and the zero-order Boltzmann equation with isotropic collisions is satisfied by any iso- 
tropic distribution function. From equation (A3) if v = 0, one can also see that the solu- 
tion then becomes ^o( v H ,v l)’ case i s n °t considered in this report. This same 

technique can be applied to the zero-order Boltzmann equation with anisotropic collisions 
(eq. (81)) to show that its solution is also any isotropic distribution function fg(|v |). 


74 


IB Ml 


■ III I I III 


I 


I I 



APPENDIX B 


DERIVATION OF FIRST-ORDER DISTRIBUTION FUNCTION 
FOR ISOTROPIC SCATTERING 

In this appendix the equations 

, , df. 

i(w u - k . v)f 1+ — = h(v) 


and 


Mii=-(7xs c ) 

dr v c/ 


are solved to obtain the solution given by equation (26). By direct substitution, 
shown that equation (B2) is satisfied by 

v(r) =Dp c (r 0 - T)J . v(t 0 ) 

where 

D i; j(<p) = n t nj + (6^ - ninj^cos cp + ejj k n k sin cp 


To solve equation (Bl), first find the solution of the homogeneous equation 

/ - X df TT 

»K- k + 

Rearranging yields 

d in f H _ . 

-= - k • v) 


dT 


and the solution of this equation is given by 


f H ( T ) = expjjia^T + i § k • v(T')drJj 


(Bl) 

(B2) 

it is easily 
(B3) 

(B4) 

(B5) 

(B6) 

(B7) 


Note that |v(r)| is not a function of t as can be shown from equation (B3), and hence 
the term a)^(v) = oj - if(v) may be extracted from the integral. Now by using the 
method of variation of parameter, it is assumed that 
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fl(f ( T )] = C(t) fjj(r) 


(B8) 


is a solution of equation (Bl) and, as a result, 




(B9) 


After integrating, 


C(T) = y 


_ r r h|f(T’g 


f H ( r ’) 


dr’ + Constant 


and equation (B8) then becomes 


(BIO) 


fl[v (r)] = <jconstant + § dr^f^r) 

Let r = Tq, solve for the constant, and obtain 

_ hM . ro dT , 


Constant 


ol _ f o h ly( r 

f ”( r o)"^ f w (r,) 


L H 


H 


(Bll) 


(B12) 


which when substituted back into equation (Bll), yields 




%< T > 


dr' 


Then using equation (B7) yields 


(B13) 


f H (T) 


= exp 


-iw^r - Tjj + i y k • v(r*)dr’ 


(B14) 


But Tq is arbitrary so let Tq be and since co is on a Laplace integral path, it 
can be seen from equation (B14) that 


f H (r > 

e h( t o) 


- 0 


(B15) 
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APPENDIX B - Continued 


Equation (B13) then becomes 


f 1 = y h[V (t')| exp j^iWj/r - r') + i ^ * IT • v(r")dr^JdT' 


Let 


t = T - T* 
t’ = T - T" 

and from equation (B3) 

v(t') = D ( w c t) * v(t) 
v(t") = D(w c t’) • v(t) 

Equation (B16) then becomes 

f i = y h[jD (^ct) • ^jexp -iw^t + i ^ k • D(o> c t') • v dt' 

The integral in the exponential is evaluated and a new tensor is obtained 

pt 

Lij^ct) h ca c ^ Dij (a> c f) dt’ = n i n j w c t + (6^ - n^sin u> c t - e ijk n k (cos w c t 
and then by defining cp = co c r, equation (B19) becomes 
f-^(v) = e“^^’ v ^h[p(<^) • \Tjd<p 

where 


dt 


co cp - k • L • v" 

*(<P,v) = i — 




From equation (21), 


h(v) = £ Ej • ^ f f x (v ’ ')dO* 

m 1 dv 4n J 1 


(B 16 ) 


(B17) 


(B18) 


(B19) 


l) (B20) 


(B21) 


(B22) 
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APPENDIX B - Concluded 


III 


and h[p(<p) • \^J must now be found. The integral term is a function of |v| 
so is unchanged. Also, f^ is a function of |v"|; thus, 

9f 0 (v > _ v 9f 0 _ v f . 

3v v av v 0 

Equation (B22) then becomes 

h[D(<rt • v] = ^ Sj • D • v + £ J f,(v ')dQ' 
and substituting equation (B24)into equation (B21) yields 


f' 

1 n 


f i (v) = — (* e-*(<P>v) ± g • d • v d<p + -0— f e-^(^ v )d<p f f-OT'JdO 

1 w c J 0 ml v 4nw c J 0 J 1 

Taking the integral ^ dfi' of this equation and rearranging yields 

r e [ A(v’)dO’ 

J f 1 (v , )df2’ = ^ 


mw, 


1 - 


47TUJ J ^0 


yy e -*(<P,v') d n' 


where 



D • v ’ 1 J. dtp 


Substituting this equation back into equation (B25) leads to 


fj(v) 


e 

mw c 


A(v) + 


4 ^ [ A < r ’> dn ’ £_ ** 

1 f fe'Wd?) dJ2* 
47TW C j Jq 


Thus, equation (26) has been obtained as desired. 


only and 


(B23) 


(B24) 


(B25) 


(B26) 


(B27) 
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APPENDIX C 


MAXWELL’S EQUATIONS AND WAVE DISPERSION RELATIONS 
Maxwell’s equations may be written 


e Q V • E = p 

(Cl) 

V • B = 0 

(C2) 

V X E = 

(C3) 

at 

1 _ ^ 5 -r aE 

-VXB = , + « 0 - 

(C4) 


The Fourier and Laplace transforms of these equations are taken as defined earlier and 
the following equations are obtained: 


-ik • ¥. = — (C5) 

e 0 

k • B = 0 (C6) 

k XE = wB (C7) 

_i ^f^ = T + iwge o (C8) 


Electrostatic Approximation 

First, consider what is commonly called the electrostatic approximation. This 
approximation is valid for waves whose electric field is very nearly parallel to the prop- 
agation vector k. The dispersion relation can now be derived from equation (C5). The 
charge density in the plasma is given by 

P = Pi + P e - en J fj(v)d 3 v (C9) 

However, for the case being considered, the ion density p^ is exactly canceled by the 
zero-order electron density p Q . The part that remains is a polarization charge density 
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that can be written as the gradient of a polarization vector 


P = 


-V • 



(CIO) 


Now define a polarization tensor by P = v • E. Taking the transforms as defined earlier 
yields 


ik • 7r • E = p = 



(Cll) 


and from this equation u may be obtained when given f j. Substitute equation (Cll) into 
equation (C5) to obtain the dispersion relation 


k • + i-j . E = 0 (C12) 

where I is the unit tensor. It was assumed earlier that k and E were parallel and 
at this point it is further assumed, without loss of generality, that k = ke^. Then, 

Jr 11 

D(o»,k) = 1 + -ii = 0 (C13) 

0 

This equation is the dispersion relation used to describe electrostatic or longitudinal 
waves. In a plasma Landau waves (without magnetic field) and cyclotron waves (with 
magnetic field) are in this class. 


Electromagnetic Waves 

For this case nothing is assumed about the E and B vectors of the waves. How- 
ever, the perturbation distribution function fj(v) is still assumed to be known, and by 
using 


J = -en y vf^vjd^v (C14) 

an expression for the current density in the plasma is obtained. Then the conductivity 
tensor a defined by 

J = <x ■ E = -en ^ vf^(v)d^v (C15) 

is obtained. Next combine equations (C7) and (C8) to obtain the equation 
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i + r-2— + (-) 2 (Cir - k 2 f)| • E = 0 

io>e 0 \u)J \ ' 

But for this to be true, the equation 

D (a,,u) = |i + 1^- + (s.) 2 (k r - k 2 r^ 


= 0 


(C16) 


(C17) 


must first be satisfied. Equation (C17) gives the relationship between oj and k and 
is the wave dispersion relation describing the waves propagating in the medium whose 
conductivity tensor is a. 



APPENDIX D 


SIMPLIFICATION OF CONDUCTIVITY TENSOR FOR k 1 B 


In this appendix the general conductivity tensor given by. equation (34) is reduced to 
the form given by equation (42). This result is exact and no approximations are made. 
The steps needed to obtain equation (42) necessarily involve large amounts of algebraic 
manipulations of a generally straightforward nature; thus only a brief outline of the pro- 
cedure is given here. First, G(v) given by equation (36) is evaluated, and from 
appendix E 


G(v) = — 

Cx> 

V 


1 + S 


2 1 + 1 


Equation (34) then becomes 


(Dl) 


a ij 


2 oo -i 

^ 7^0 Vf 0 ( v) I Vi A.(v)dfi dv 


U) f _ r>°° 

- — 2 y ™/|)W 

47TCt>2 v U 


y Aj(v)dfi y y v^ ^d^ d^ 


a> + i^S 


2Z + 1 


dv 


(D2) 


Now consider each component of the tensor in turn. First cr^g = = Cg^ = Ogg = 0. 

Next consider crgg where it can be seen that the second term vanishes, and by using an 
integral given in appendix E 


2 P°° y3f o 11 

<r„ = 47riw“e n \ -/± + S 

33 P 0 J n & )2 


Rdv 


(21 + 1 )( 2 1 + 3)JJ 


(D3) 


is obtained. For and °21 again integrals worked out in appendix E are used to 

find 


4 ™ e o w ? w c r vf o s (2FTT 


ct 12 " " ct 21 


I 


® cu + ivsf ^ 


dv 


21 + 1 


(D4) 
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r - 


For a 22 , 



Finally, for a jj, 



but after rearranging one can obtain 

iwe n 47ra>^ r»°° 


CT 11 ~ 


: q^^p r 

k 2 ^0 


31 + 1 


a) + i^S 


21 


±4 dv 


(D5) 


(D6) 


(D7) 


The algebra necessary to derive these equations is quite extensive but generally straight- 
forward. One can show the equality of equations (D6) and (D7) most easily by subtracting 
the two and showing that zero is obtained. At one point it is helpful to make use of the 
identity 


D • (IT • L) s L • k (D8) 

This equation can easily be proved by using equations (29) and (30) for the tensors and 
showing that both sides yield the same result. 
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APPENDIX E 

THE EVALUATION OF SEVERAL PERTINENT INTEGRALS 

In this appendix several integrals used in appendixes D and F are evaluated and the 
rest are stated without proof as they are all done in a similar manner. 

Consider first the integral 

Ij = e' #(<p ’ v ^dfi (El) 


where 

M<P,v) = irjcp - - 

w c 

and 

V = — 

(jl) 

C 

Let F = ^ — and rewrite equation (El) to obtain 

Ij = e~ i? ^ y e 1 ^ ’ ^dfi (E2) 

This integral is easily evaluated and yields 

Ij = dS2 =-e -17?,?? 47rj 0 (Fv) (E3) 

where Jq(Fv) is a spherical Bessel function of the first kind given in terms of the ordi- 
nary Bessel function of the first kind by 

] n (z > + 

Next consider the integral 

! 2 = J v e_$ ^’^)dO (E4) 


dak 
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which can be written as 


I = _i _L C e-*^) dS 2 = _i 4_ i 

* 9F J 8F 1 

Then the following equation is easily obtained: 

Ig = y ve‘^^ ,v ^ dfi = 477vie"^ <;0 j j(Fv) ^ 


(E5) 


(E6) 


In a similar manner, it can be shown that 


I 


v.v,e“^ > ^» v ) dfi = 4ire - * 7 ^ 
J ^ 


vji(Fv) F.F 


I,i F 


jj- T V Fv > 


Now investigate integrals of the type 


fe'^FWd? 

J 0 


(E7) 


(E8) 


In all the integrals of this type needed in appendixes D and F, it is possible to expand the 
F(<p) in powers of sin (<p/ 2) and then to evaluate the integral in terms of a series. To 
dembnstrate, consider 


I 3 = y e-*^) d(p 

After using equation (E3), equation (E9) becomes 
I 3 = § e _1, 7^ 477j 0 (Fv)d<p 


(E9) 


(E10) 


But for the case being considered, that is, lT = ke^, one can show F = (2k/o> c )sin(<p/2). 
After expanding ]q(Fv) in its series representation, equation (E10) becomes 


I = 4 tt f (-/ C°° e^sin 21 % d<p) 

3 ^ (21 + i):\ w c j \ J 0 2 / 


(Ell) 


The integrals in equation (Ell) are Laplace transforms and were evaluated by Bateman 
(ref. 52). Using his results yields 
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I 3 = f C e-*(<P& d <p dO = + 

3 J J 0 iu v \2l + l) 


where, for convenience, the symbolic notation 


■ih) - 1 


( kv / W c) 


2 1 


tl (^ 2 - 1 ) ...U 2 -* 2 ) 


(E12) 


(E13) 


has been used where 


„ = to - Mv) 
w c 

The remaining integrals are solved in a similar manner; therefore, they will be given 
here without proof. The list of integrals follows: 


y e"^dO = 4rj Q (Fv)e' 


ir](p 


y ve'^dfi = 4?rivj j(Fv)e -i? ^ <;f ’ ^ 


y 


v-v^e-^ dfi = 4 it 

J l 


J,(Ft) ff 

\ ? -v- - -y v j 2 (ft) 


CC e- # d^ dfi = 1 + s(-J— ) 

J Jo iu > u L \2l + l[ 

y 0 e - iW c° S f 

■ „ M Fv ) 


-irj(p 


r e -iw!^ # .i { tS r. 1 1 

Jo F i7?\3 |_(2i + 1)(2Z + 3)_ 


(E14) 
(El 5) 

(E16) 

(E17) 

(E18) 

(E19) 

(E20) 


I 0 e-iWj^Fvtd^.-is- 


21 


(2 1 + 1) (22 + 3) 


(E21) 
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e~ 171<p sin 2 ^ j 2 (Fv)d<? = ^ 



j* e~ lv<p sin ^ cos ^ j 2 (Fv)d<p = 




(E22) 

(E23) 

(E24) 

(E25) 






APPENDIX F 


COLLISIONAL CYCLOTRON DISPERSION RELATION FROM 
ELECTROSTATIC APPROXIMATION 

In appendix C it is shown that for waves that are nearly longitudinal, the electro- 
static approximation yields the dispersion relation 

D(co,k) = 1 + — = 0 (FI) 

e 0 


where 


p = -en y f j(v)d^v = ik • rf • E (F2) 

From equation (26), an equation for the first-order distribution function fj is obtained 
which when used in equation (F2) yields 



(F3) 


Recalling that B = Beg and it = ke^ and using the integrals given in appendix E 
results in 


^1 = ^0 


m r 

\k / J 0 


'°° vfQG(v)(o> - ip) 


u> + iuG(v) 


dv 


(F4) 


where 



(*v/«c)” 

(21 + l)(v 2 - l) . . . (t] 2 -l 2 ) 


and 




88 



APPENDIX F - Concluded 


Substituting equation (F4) into equation (FI) yields 


w \ vf 0 G(v)co 

D(u> k) = 1 + 47r(-j-£) \ . . . 

\ k / Jq w + u/G(v) 


dv 


(F5) 


and the dispersion relation given in equation (51) has been found but in a much simpler 
manner. The more complicated derivation was used earlier to show how the electro- 
static dispersion relation could be obtained from the more complete theory. 
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QUADRATURE RULE FOR WEIGHT FUNCTION exp(-x 2 ) ON 


In several of the dispersion relations solved in this report, it was necessary to 
evaluate integrals numerically. These integrals were all of the form 

f e" x2 f(x)dx (Gl) 

J 0 

and were performed by a Gauss quadrature rule. This rule enables us to write 

oo N 

1 ZVft) < G2 > 

0 1=1 

It is easily shown that if w^ and x^ have been chosen correctly, this equation is exact 
if f(x) is a polynomial of order less than 2N. From Gallant (ref. 53) and Wilensky 
(ref. 18) the weights w^ and the nodes x^ are used for N = 20 and 40. These values 
are given in table II. 

TABLE II.- WEIGHTS AND NODES FOR THE QUADRATURE RULE USING THE WEIGHT 

FUNCTION -p exp(-x 2 ) ON [o,o0 
V K 

(a) 20-point rule 


x 

l 

o. 1427950969991826E-01 
0. 7U63130039219270E-01 
0.1808615630580397E 00 
O.32943335606J+2905E 00 
0.5160505430615322E 00 
0. 736255457580891 4E 00 
0.9858735750375285E 00 
0.126128 90161027 67E 01 

0.1559579645209661E 01 

0.187856191930298IE 01 
0. 2216794165387 664e 01 
0.2573577820826277E 01 
O.2948989746787235E 01 
0. 3343981379861640E 01 

0.3760599932910770E 01 
0.4202442600421132E 01 
0. 4675608847794480E 01 
0.5190901686974998E 01 
0. 576998 5165567762E 01 
0.647055838706457^E 01 


0. 4126248095354063E-01 
0. 9401377007261654E-01 
0.1398673127403206E 00 

0.1704175033159441E 00 
0. 17649482600028 38E 00 
0.1545791839772282E 00 

0.1123067291104619E 00 
0. 6604828742706440E-01 
0. 3061944925866195E-OI 

0. 1088291653324355E-01 

0. 2879302149689494E-02 
0. 5488989148150679E-03 
0. 7261283758335416E-04 
0. 6366757009378284E-05 
0. 348763554069259^E-06 
0. 1100893097830079E-07 
0. 1778423417441857E-09 
0. 1214074765543346E-11 
0. 2448428821057576E-14 
0. 59930 73454594019E-18 
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TABLE II.- WEIGHTS AND NODES FOR THE QUADRATURE RULE USING THE WEIGHT 
FUNCTION JL exp(-x^) ON (o,J] - Concluded. 

(b) 40-point rule 


0. 51500103 5639331^E -02 
0. 270T80249l*829ll3E-01 
0. 6629858021142996E-01 
0. 122439601 5988346E 00 

0.1949786809306882E 00 

0. 2832755895057396E 00 
0.3866018422176653E 00 
0. 504170513 4328489E 00 
0.6351645938060005E 00 

0.7787623792828677E 00 
0.9341589652282645E 00 

0.1100583488138671E 01 
0.1277312213015507E 01 
0.1463677876957118E 01 
0.1659075877002548E 01 
0.1862967959865027E 01 
0.2074884066030345E 01 
0. 229442293170663 5E 01 
0.252125198 4268 593E 01 
0.2755106998376953E 01 
0. 2995791923 381212E 01 
0.3243179256519709E 01 
0.3497211327571455E 01 

0.3757902886114078E 01 
O.4O25345452403971E 01 
0. 42997140233 22116E 01 

0. 458127 6943399332E 01 

0. 4870410105599751E 01 
0.5167617222469367E 01 

0. 5 1 )-7355886ll03993E 01 
O.578909456041786OE 01 

0.6115345232898952E 01 
0.6453788 427 291654E 01 
0.6806409661424949E 01 
O.7175955680231676E 01 
0.7566388429073037E 01 
0.7983778287554144E 01 
0.8438 3115413 5517 3E 01 
0 . 89498 413 55010714E 01 
0.9570481178890405E 01 


0. 1490554725H8071E-01 
0. 3452733473086753E-01 
0. 5365120033130722E-01 

0. 71628 48 48276338 9E-01 

0.8748201803249430E-01 
0. 9990783959338579E-01 

0.1074557174139692E 00 
0.1088800688850989E 00 

0.1035938920774140E 00 
0. 9206274103552843E-01 
0. 7592560230950505E-01 
0. 5769699H393027OE-OI 

0. 4009918101795804E-01 

0. 2529440279616882E-01 
0. 1437098784698673E-01 
0. 729769626394253IE-02 
0. 3286876449615393E-02 
0. 1302941977664821E-02 
0. 4510407826591482E-03 
0. 1352649280556334E-03 
0. 3485330737193931E-04 

0. 7649469079803097E-05 
0. 141694OO72635O54E-05 
0. 2193302752467988E-06 
0. 2806496911831388E-07 

0. 293309698298626IE-O8 

0. 2470041550273338E-09 
0. I650361609485315E-IO 
0.8593168390112723E-12 
0. 3413713673581817E-13 
0. 1008670314083280E-14 
0. 2148684677740992E-16 
0. 31731788 11113663E-18 
0. 3088596815526213E-20 

0. 1851864065560540E-22 

0. 6221943052492195E-25 
0. IOI7307376257966E-27 
0. 643044527 9019888E-31 
0.1017 466944507 922E - 3 4 
0. 134940456968 5499E-39 


APPENDIX H 


ASYMPTOTIC EXPANSION OF G(v) 


It has been pointed out that one can write 
G(v) = -1 + 2 F 3 

Meijer (ref. 31) gives the theory of the asymptotic expansion of the generalized hyper- 
geometric series p F^, and after specializing his results to the 2 F 3 in equation (HI), 
one obtains 


1, 1/2; 3/2, !-??,! + 




(HI) 


G(v) = 


jvfif 


2p sin rjTt 


A cos 



+ B sin 



+ 



(2 1 + l)p 5 


1 


(H2) 


where 



_ co - iv(v) 



The h terms are given by the recurrence relation 
1 ~ a l-l 1 + 1 


(H3) 

(H4) 


(H5) 


(H6) 


(H7) 
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where 


hj = 1 (H8) 

v ; - ^ + B ;-i(l + !) ♦ B i+i(| + 1) < H9 > 



(l odd) 


(l even) 


and B^(x) are the Bernoulli polynomials. 


(H10) 
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NUMERICAL CONTOUR INTEGRATION 


Let us consider the contour integral of some function F 


= y F(z)e“ 


z2 dz 


(ID 


with the contour given by 


z = x 1 + 


2iyx m 

x 2 + x 2 ‘ 
x m + x 


( 12 ) 


where 


y = 


T m 




0 ^ x < 00 


(13) 


Note that x is real, and one can easily show that Max[lm(z)J = y m at x = x m . A 
plot of this contour is shown in figure 28. Substituting equation (12) into equation (II) 
yields 


= £ F(z)e- 


z 2 dz 


-f' 


x 1 + 


2i yx 


m 


x^ + x 2 


_Y" 

e exp 


4yx 2 x^ 


yx 


x2 + x m 


-1 + 


m 


X 2 + X 2 

m> 


. 2irx m( x m - x2 j 

-*■ / ft n v ft 


(4 - * 2 ) 2 


dx (14) 


From appendix G, the Gauss quadrature rule is 

N 


f f(x)e- x2 dx = & £ Wjf(xj) 
U 1=1 


(15) 


94 



APPENDIX I - Concluded 



Figure 28.- Normalized mapping of contour in complex z plane. 

Note that Re(z) = x. 


and applying this rule to equation (14) yields 

N 


f F(z)e- 2 cfc = & F(z,) 


1 = 1 


where 


2iyx: 


m 


z ' = Xill + 4^F 


™i = w i 


’ , “<(4 - 4 ) 

j A n v n 


( 


x m + x ?) 2 


exp 


4yx?x^ 
l m 


X m + 

m L 


-l + 


yx m 

x m + x ? 

m L 


(16) 


(17) 


(18) 


By using the technique described, it is possible to do very accurate numerical contour 
integration by means of the Gauss quadrature rule. Note that no rapidly oscillating 
terms are introduced in the integrand of equation (14). For best results it was found 
that y should be kept small, for example, y < 0.5. 


APPENDDC J 


DERIVATION OF TRANSFER COLLISION CROSS SECTION Q^(v) 
From the definition. 




Y lm (e,<p) 


<W V > = y 0,v >*> 

and use of the Legendre polynomial expansion 

OO 

o-(v,y) = ^ cr n (v) P n (cos x) 


dO' 


n=0 


yields 


Q fm = Z a n (v) ^ P n (cosx) 


1 - 






dn’ 


From figure 29 it can be seen that 


cos x = cos 9 ' cos 8 + sin 6' sin 9 cos {cp' - (p) 


which allows the use of the addition formula 

n 


5 n (cos x)= l (-) k P^(cos 9') P” k (cos 0) e im ^ T ' 


<P) 


k=-n 


Then evaluate the integrals in equation (J3) and obtain 


< %m = 4ff l <T n - 


0 21 + 1 


(Jl) 


(J2) 


(J3) 


(J4) 


(J5) 


(J6) 


This equation is identical to that obtained from substituting equation (J2) into Allis' 
transfer collision cross section 


Qj(v) = 2vr j* a(x,v) [l - P^(cos x)j sin x dx 


(J7) 
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Figure 2 9-- Geometry of collision velocities v and v' before 
and after impact, respectively. 


Thus it is concluded that 

<W v) = Q Z (v) 


(J8) 


for all m as claimed in the text 



APPENDIX K 


CONDUCTIVITY TENSOR FOR k || B WITH 
ISOTROPIC SCATTERING 


In this appendix the conductivity tensor from equation (34) will be obtained for the 
case of propagation parallel to the magnetic field. By taking k || B, that is, k = ke^, 
equation (38) yields 

co - kvo 

H<P,v) = i -^7 ^ <P (Kl) 

w c 

Substituting equation (Kl) into equation (34) yields 
a 13 = a 23 = a 31 = °32 = 0 


a 


a 


a 


Vo r £vr 

co c j V \J 0 


ii " °22 zr. 


e~® cos cp d<p) d 3 v 


12 - - a 2 i = ^7 1 nr(f 0 e ~* sin f d<p ) rv 




33 


c ^0 


r S t dn ]' 


47TCO, 


i . _ii_ c r d<p do* 

4ttco c J Jq ^ 


' dv 


(K3) 

(K4) 


(K5) 


Upon evaluating the integrals in equations (K3) to (K5), one finds 


2 

f«vr 


2 r I o v i 
a ll " a 22 " Vo J T~ 


- kv 3 


"1 - K - kv 3)' 


d 3 v 


(K6) 
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f'v* 


ff 12 = - ct 21 = V 


f 0 1 

0 J v 


u> 


°33 “ "■ 


“l - K - kv s ) 2 

1-^oth-U dy 

k 2 j 0 ° 1 + iii coth -i x 

kv 


d 3 v 


47rio) 3 e n w /*» 


(K7) 


(K8) 


where 



The angular part of the integrals in equations (K6) and (K7) can be easily evaluated but 
will not be needed at this time. 


APPENDIX L 


PROOF THAT THE ANISOTROPIC DISPERSION RELATION YIELDS 
THE ISOTROPIC DISPERSION RELATION IN THE LIMIT 
OF ISOTROPIC SCATTERING 


In this appendix the dispersion relation expressed by equation (108) is reduced to 
the one in equation (116) in the limit of isotropic collisions. For isotropic collisions, 
that is, a = cr(v), it can be seen from equation (88) that all the terms are equal. 
Hence let = v(v) for all l and equation (108) then yields 

/ , \2 2 Trct)^ n°° 

= 1 - (f) + -7ST J 0 g i (v> v2f o (v) dv (L1) 

where 

2 3 8 n 2 - 1 

g ± 3\ ± - 5X ± - 7\ ± - ' ' ‘ - (2n + 1)\ ± 

and 




U) ± (j0 

V c 
kv 


Take now the part of the integrand in brackets in equation (116) and expand it in a series 
in powers of X" 1 to obtain 




(> - x D 


coth“l X = -- 

* V 


>-2 \ -4 

1 X ± ± 

'3 + 3T + 5 T + 


-2n 


• • + 


(2n + l)(2n + 3) 


(L2) 


Convert this asymptotic series into a continued fraction by the method of Viskovatoff as 
described by Khovanskii (ref. 43) to obtain 




coth - -*- \ . 


2 3 8 

3\ ± - 5A. ± - - 


n 2 - 1 

- (2n + 1)\ ± - ' * ‘ 


(L3) 


Hence equation (116) yields an equation identical with equation (LI). It has therefore 
been shown that equations (108) and (116) are identical in the limit of isotropic collisions. 


APPENDIX M 


SIMPLE DERIVATION OF DISPERSION RELATION 
FOR WAVES WITH k || B 


In this appendix the dispersion relation for propagation parallel to the magnetic 
field is derived in a simpler but less rigorous manner (ref. 45) than was used earlier. 
The assumption is made that the waves are transverse in the course of the derivation; 
therefore, only the dispersion relation for the transverse electric waves is obtained. 

For isotropic collisions in a Lorentz gas, the Boltzmann equation is written as 

+ v • — - — (e + f xb) - + f(v’)dS2’ (Ml) 

at 9 r m v av 4 ttJ 


The ansatz is made 

~\ 

f(v) = f 0 (v) + fj(v) 

E = 0 + E 1 > 

§ = V § 1 J 


(M2) 


and the zero- and first-order equations for fg and f^ are obtained. The zero-order 
equation has as its solution any isotropic fg(|v I) (appendix A), and the first-order 
equation is solved by means of Fourier and Laplace transforms in space and time based 
on the wave representation exp(io>t - iic • r). Taking Bq = BgSg and writing the first- 
order equation in spherical coordinates then yields 

i (u> - kv cos 0)1, + oj — - = Cl^vOdST + — E 1 • — - (M3) 

V u / 1 C 47T J 1 m J- 0y 


Expanding the dot product in component form results in 


af 


E 


0 „ l 0 


1 av 


sin 6 
~2 




(M4) 


where E ± = E x ± iEy are the fields of the right- and left-hand polarized waves, respec- 
tively. Now the solution 
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is assumed and from equation (M3) 


f, = - 


le 


E ± fg sin 9 


± 2m o> v t a> c - kv cos 6 


(M6) 


is easily obtained. From Maxwell's equations of electrodynamics, the current density is 

>2 


T= -e n y vfid 3 v= g) - 1 


icue 0 E l 


which yields 


j ± = -2 t r en ^ sin^0f ± d0 dv = (^fj^ - 1 


icue 0 E ± 


(M7) 


(M8) 


Then, substituting equation (M6) into equation (M8) yields 


Trine^E , p°° pir sin 2 0fl d0 dv 


m 


a 


0 J 0 T w c “ kv cos 0 


/ck\ 2 _ x 

\ CO / 


iwe o E ± 


(M9) 


Equation (M9) then gives the dispersion relation 




2 7 rw* r>°° pi t V 3 sin^0f 


0 J 0 Wj, ± co c - kv cos 9 


d9 dv = 0 


(M10) 


where 


w n = — 
P me. 


Evaluating the integral over 9 then yields 

D ± ("' k > = 1 - (t)‘ 


/ 2 277 CO 2 

/ 

n\ - 

( ck ) + , p \ v 2f' 

X^ + 

1 - X" coth" 1 ^^ 

V co / cok Jq 0 

± V 

±/ ± 


dv 


(Mil) 


where 


X 


± 


CO + CO 


kv 
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which is the same as equation (116). This derivation demonstrates that the integral term 
in the Boltzmann equation (eq. (79)) does not contribute to the dispersion relation for 
transverse waves propagating parallel to the magnetic field. Note the integral vanishes 
when obtaining equation (M6). From equation (M10) it can also be seen how D_(a>,k) 
and D + (ci>,k) are associated with right- and left-hand polarized waves, respectively. 


APPENDIX N 


REDUCTION OF EQUATION (116) TO FRIED FUNCTION 
REPRESENTATION IN THE APPROPRIATE LIMIT 

In this appendix the Fried function representation of the dispersion relation 
(eq. (119)) is obtained from equation (116) under the assumption of (1) a constant collision 
frequency and (2) a Maxwellian electron distribution. Integrating equation (116) by parts 
yields 


.2 4t no* 

D ± (u>,k) = !-(-) - 1 vf 0 COth dv 


If fp is Maxwellian, 


,rf - 0 f' 

" ~Y f o 


and the integral in equation (Nl) becomes 


Is I 0 vt o coth ~ lx ± Av = -Y ^ 0 f i coth_lx ± dv 


Doing a second integration by parts and substituting equation (75) for fp results in 


x 0 r e ~ x ^ d .. 

2 ^ 3/2 Jq x 2 _ 


where 


<j0 ± Cl) 

. _ V C 

°~ kv o 


The integrand is an even function; thus, 


1 = - 


4?r 0/ * V n x* - X 


2 _ >2 
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APPENDIX N - Concluded 


Expand equation (N5) in partial fractions to obtain 



Putting x x in the second integral yields 


(N6) 


I = 


4^ 


r 


a-* 2 

X - 


dx 


(N7) 


Since the dispersion relation is defined for k real and to on a Laplace integral path, 
equation (N7) is the proper equation for Ihi^Xq) < 0. For other values of \q, the 
analytic continuation is used. Hence define 


: ( x o) 


»i-f 




\fn -°° x ' x o 

and equation (N7) then becomes 

z ( x o) 


dx < 0^ 


I = -- 


47TV, 


0 


Substituting equation (N9) into equation (Nl) yields 


D ±K k > = i - (f )‘ 


tot 


u>kv. 


0 


to m to 
v c 

kv 


and hence gives equation (119). Fried and Conte (ref. 35) define 




dx 


JF y/H J -oo X - £ 

and by using equation (N8) one can show 

z ( x o) = z f( x o) 


(Im(4) > 0) 


(N8) 


(N9) 


(N10) 


(Nil) 


(N12) 


Equation (N10) can then be written in terms of the Fried function (plasma dispersion 
function) if desired. 
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